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ABSTRACT 

We present the Sloan Low-mass Wide Pairs of Kinematically Equivalent Stars (SLoWPoKES), a catalog of 
1342 very- wide (projected separation > 500 AU), low-mass (at least one mid-K - mid-M dwarf component) 
common proper motion pairs identified from astrometry, photometry, and proper motions in the Sloan Digital 
Sky Survey. A Monte Carlo based Galactic model is constructed to assess the probability of chance alignment 
for each pair; only pairs with a probability of chance alignment < 0.05 are included in the catalog. The 
overall fidelity of the catalog is expected to be 98.35%. The selection algorithm is purposely exclusive to 
ensure that the resulting catalog is efficient for follow-up studies of low-mass pairs. The SLoWPoKES catalog 
is the largest sample of wide, low-mass pairs to date and is intended as an ongoing community resource for 
detailed study of bona fide systems. Here we summarize the general characteristics of the SLoWPoKES sample 
and present preliminary results describing the properties of wide, low-mass pairs. While the majority of the 
identified pairs are disk dwarfs, there are 70 halo subdwarf pairs and 21 white dwarf-disk dwarf pairs, as well 
as four triples. Most SLoWPoKES pairs violate the previously defined empirical limits for maximum angular 
separation or binding energies. However, they are well within the theoretical limits and should prove very 
useful in putting firm constraints on the maximum size of binary systems and on different formation scenarios. 
We find a lower limit to the wide binary frequency for the mid-K - mid-M spectral types that constitute our 
sample to be 1.1%. This frequency decreases as a function of Galactic height, indicating a time evolution of 
the wide binary frequency. In addition, the semi-major axes of the SLoWPoKES systems exhibit a distinctly 
bimodal distribution, with a break at separations around 0. 1 pc that is also manifested in the system binding 
energy. Comparing with theoretical predictions for the disruption of binary systems with time, we conclude 
that the SLoWPoICES sample comprises two populations of wide binaries: an "old" population of tightly bound 
systems, and a "young" population of weakly bound systems that will not survive more than a few Gyr. The 
SLoWPoKES catalog and future ancillary data are publicly available on the world wide web for utilization by 
the astronomy community. 

Subject headings: binaries: general — binaries: visual — stars: late-type — stars: low-mass, brown dwarfs — 
stars: statistics — subdwarfs — white dwarfs 



1. EvlTRODUCTION 

The formation and evolution of binary stars remains one of 
the key unanswered questions in stellar astronomy. As most 
stars are thought to form in multiple systems, and with the 
possibility that binaries may host exoplanet systems, these 
questions are of even more importance. While accurate mea- 
surements of the fundamental properties of b inary systems 
provi de constraints on evolutionary models (e.g. lStassun et alJ 
l2007h . knowing the binary frequency, as well as the distribu- 
tion of the periods, separations, mass ratios, and eccentrici- 
ties of a large ensemble of bi nary systems are crit ical to un- 
derstanding binary formation ( iGoodwin et alj2007l and refer- 
ences therein). To date, multiplicity has been most extensively 
studied for the relativ ely bright high- and solar-m ass local 
field p opulations (e.g. iDuquennoy & Mayoilll991l hereafter 
IDM91I) . Similar studies of low-mass M and L dwarfs have 
been limited by the lack of statistically significant samples 
due to their intrinsic faintness. However, M dwarfs consti- 
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2013 hereafterlBlob' and significant ly influe n ce its properties 
Since the pioneering study of iHeintzl (119691) . binarity 
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for brown dwarfs (BDs; iBouy et al.l [20031; i Close et al.l 12^531; 
iGizis et al.1 [20031; iMartm et alj 120031) . This decrease in bi- 
narity with mass is probably a result of preferential de- 
struction of lower binding energy systems over time by dy- 
namical interactions with other stars and molecular clouds, 
rather than a true represen tation of the multiplicity at birth 
( IGoodwin & Kroupal 120051) . In addition to having a smaller 
total mass, lower-mass stars have longer main-sequence (MS) 
lifetimes (Laughlin, Bodenheimer, & Adams 19971) and, as 
an ensemble, have lived longer and been more affected by 
dynamical interactions. Hence, they are more susceptible 
to disruption over their lifetime. Studies of young stellar 
populations (e.g. in Taurus, Ophiucus, Chameleon) appear 
to support this argume nt, as their multipUcity is twice as 
high as that in the field dLeinert et al.lll993l;lGhez et al.ll 19971; 
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iKohler & Leinerllll998l) . However, in denser star-forming re- 
gions in the Orion Nebula Cluster and IC 348, where more 
dynamical interactio ns are expe cted, the multiplici ty is com- 
parab l e to the field (ISimon. Close & Be ck 199§, iPetr et al.l 
Il998t iDuchene. Bouvier. & Sirnonill999h . Hence, preferen- 
tial destruction is likely to play an important role in the evo- 
lu tion of binary systems. 

IDM9 ll found that the physical separation of the binaries 
could be described by a log-normal distribution, with the peak 
at a ^30 AU and crioga ^1-5 for F and G dwarfs in the lo- 
cal neighborhood. The M dwarfs in the local 20-pc sample 
of lFM92l seem to follow a similar distribution with a peak at 
a ^^3-30 AU, a result severely limited by the small number of 
binaries in the sample. 

Importantly, both of these results suggest the existence 
of very wide systems, separated in some cases by more 
than a parsec. Among the nearby ( d < 100 pc) solar- 
type s tars in the Hipparcos catalog, iLepine & Bongiomol 
(120071) found that 9.5% have companions with projected or- 
bital separations s > 1000 AU. However, we do not have 
a firm handle on the widest binary that can be formed 
or on how they are affected by localized Galactic poten- 
tials as they traverse the Galaxy. Hence, a sample of 
wide binaries, especially one that spans a large range of 
heliocentric distances, would help in (i) putting empirical 
constraints on the widest binary systern s in the field (e.g. 
Reid etal.l j2001at |Burgasser et al.l |2003l 120071: IClose et al.l 
2003, 2007|), (ii) understanding t he evolution of wide bina- 
ries over time (e .g., Wei nberg. Shapiro. & W asserman 1987; 
iJiang & Tremaine2009) . and (iii) t racing the inhomogeneities 
in the Galactic potential (e.g. iBahcall. Hut. & Tremaind 
19851 iWeinberg et al.l [19871 lYoo. Chaname. & Gouldl 120041 
Ouinn et al.ll2009l) . 

Recent large scale surveys, such as the Sloan Digital Sky 
Survey (SDSS; York et al. 2 000), the Two Micron All Sky 
Survey (2MASS:ICutri'eran i2(ffl and the UK IRT Infrared 
Deep Sky Survey (UKIDSST lLawrence et akl 120071) . have 
yielded samples of unprecedented numbers of low-mass stars. 
SDSS alone has a photo metric catalog of more than 30 million 
low-mass dwarfs (iBlOt) . defined as mid-K - late-M dwarfs for 
the rest of the pa per and a spectroscopic catalog of more than 
44000 M dwarfs (IWest et al..2008.) . The large astrometric and 
photometric catalogs of low-mass stars afford us the opportu- 
nity to explore anew the binary properties of the most numer- 
ous constituents of Milky Way, particularly at the very widest 
binary separations. 

The orbital periods of very wide binaries (orbital separation 
a > 100 AU) are much longer than the human timescale {P = 
1000 yr for Mtot= 1 Mq and a = 100 AU). Thus, these sys- 
tems can only be identified astrometrically, accompanied by 
proper motion or radial velocity matching. These also remain 
some of the most under-explored low-mass systems. With- 
out the benefit of retracing the binary orbit, two methods have 
been historically used to identify very wide pairs: 

(i) lBahcalT& Soneiral (119811) used the two-point correla- 
tion method to argue that the excess of pairs found 
at small separations is a signature of physically as- 
sociated pairs; binarity of some of these systems 
was later confirmed by radi al velocity observations 
(iLathamet all 11984. See iGarnavichI (Il988h and 
IWasserman & WeinbergI (119911) for other studies that 
use this method. 

(ii) To reduce the number of false positives inherent in 



the above, one can use additional information such as 
proper motions. Orbital motions for wide systems are 
small; hence, the space velocities of a gravitationally 
bound pair should be the same, within some uncer- 
tainty. In the absence of radial velocities, which are 
very hard to obtain for a very large number of field stars, 
proper motion alone can be used to identify binary sys- 
tems; the resulting pairs are known as common proper- 
motion (CPM) doubles. iLuvtenl (119791 [1988) pioneered 
this technique in his surveys of Schmidt telescope 
plates using a blink microscope and detected more than 
6000 wide CPM doubles with ^ >100 mas yr~^ over 
almost fifty years. This method has sinc e been used 
to find CPM doubles in the AGK 3 stars bv lHalbwachsl 
dl986'), in the revised New Luyten T wo-Tenths (rNUTT; 
[Salim & Gouldl 120031) catalog by IChaname & Gouldl 
(P004), and among the Hipparcos star s in the Lepine- 
Shara Proper Moti on-North (LSPM-N;|L epine & Sharal 
120051) catalog by iLepine & Bongiorno (2007). All 
of these studies use magnitude-limited high proper- 
motion catalogs and, thus, select mostly nearby stars. 

Mor e recently, ISesar. Ivezic. & Jurid (120081 hereafter 
'SIJ08 !) searched the SDSS Data Release Six (DR6; 
lAdelman-McCarthv etani2008l) for CPM binaries with angu- 
lar separations up to 30" using a novel statistical technique 
that minimizes the difference between the distance moduli 
obtained from photometric parallax relations for candidate 
pairs. They matched proper motion components to within 
5 mas yr^^ and identified ^22000 total candidates with ex- 
cellent completeness, but with a one-third of them expected 
to be false positives. They searched the SDSS DR6 cat- 
alog for pairs at all mass ranges and find pairs separated 
by 200 0-47000 AU, at distances up to 4 kpc. Similarly, 
iLonghi tano & Binggeli ( 2010) used the angular two-point 
correlation function to do a purely statistical study of wide 
binaries in the -^675 square degrees centered at the North 
Galactic Pole using the DR6 stellar catalog and predicted that 
there are more than 800 binaries with physical separations 
larger than 0.1 pc but small er than 0.8 pc. As evidenced by 
the large false positive rate in lSIJ08L such large-scale searches 
for wide binaries generally involve a trade-off between com- 
pleteness on the one hand and fidelity on the other, as they 
depend on statistical arguments for identification. 

Complementing this type of ensemble approach, a high- 
fidelity approach may suffer from incompleteness and/or 
biases; however, there are a number of advantages to a "pure" 
sample of bona fide wide binaries such as tha t presented in 
this work. For example, iFahertv et al.l ( 120101) searched for 
CPM c ompanions around the brown dwarfs in the BDKP 
catalog ( IFahertv et al.ll2009l) and found nine nearby pairs; all 
of their pairs were followed up spectroscopically and, hence, 
have a much higher probability of being real. As mass, 
age, and metallicity can all cause variations in the observed 
physical properties, e.g. in radius or in magnetic activity, their 
effects can be very hard to disentangle in a study of single 
stars. Components of multiple systems are expected to have 
been formed of the same material at the same time, w ithin a 
few h und red thousand years of each other (e.g. White & GhezI 
20011: [Goodwin Whitworth. & Ward-Thompsod 12004 
Stassun et alj |2008[k Hence, binaries are perfect tools 
for separating the effects of mass, age, and metallicity 
from each other as well as for constraining theoretical 
models of stellar evolution. Some examples include 
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bench m arking stellar evolutionary t r acks ( e.g . I White et al 
19991: IStassun. Mathieu. & Valentil 120071: IStassun et al 



20081) , investigating the a ge-activity relati ons of M dwarfs 
(e.g. ISilvestri Hawley. & Oswalt 120051) . defining the 
dwarf- subdwarf bo undary for spectral classification (e.g. 
iLepine et al.l 120071). and caUbrating the metall icity indices 
jWoolf & W allersteini l200l iBonfils et al.ll2005l) . Moreover, 
equal-mass multiples can be selected to provide identical 
twins with the same initial conditions (same mass, age, 
and metallicity) to explore the intrinsic variations of stellar 
properties. In addition, wide binaries (a > 100 AU) are 
expected to evolve independently of each ot her; even thei r 
disks are unaffected by the distant companion ( IClarkelfT992l) . 
Components of such systems are effectively two single stars 
that share their formation and evolutionary history. In essence 
they can be looked at as coeval laboratories that can be used 
to effectively test and calibrate relations measured for field 
stars. Finally, as interest has grown in detecting exoplanets 
and in characterizing the variety of stellar environments in 
which they form and evolve, a large sample of bona fide wide 
binaries could provide a rich exoplanet hunting ground for 
future missions such as SIM. 

In this paper, we present a new catalog of CPM doubles 
from SDSS, each with at least one low-mass component, iden- 
tified by matching proper motions and photometric distances. 
In §|2]we describe the origin of the input sample of low-mass 
stars; §[3] details the binary selection algorithm and the con- 
struction of a Galactic model built to assess the fidelity of each 
binary in our sample. The resulting catalog and its character- 
istics are discussed in §|4] We compare the result of our CPM 
double search with previous studies in §|5]and summarize our 
conclusions in §|6l 

2. SDSS DATA 
2.1. SDSS Sample of Low-Mass Stars 

SDSS is a comprehensive imaging and spectroscopic sur- 
vey of the northern sky using five broad optical bands (ugriz) 
between -3000-lOOOOA using a dedica t ed 2.5 m te l escop e 
(lYorket al.l llOOOl; IFukugita et all fl996l: iGunnetan [1991). 
Data Release Seven (DR7; lAbazaiian et al.l 12001 " contams 
photometry of 357 million unique objects over 11663 deg^ 
of the sky and spectra of over 1.6 million objects over 9380 
deg^ of the sky. The photometry has calibration errors of 
2% in u and '^1% in griz, with completeness limits of 95% 
down to magnitudes 22.0, 22.2, 22.2, 21.3, 20.5 and satura- 
tion at magnitudes 12.0, 14.1, 14.1, 13.8, 12.3, respectively 
jGunn et al. 1998|. We restricted our sample to r < 20 where 
t he S DSS/USNO-B proper motions are more reliable (see 
§ 12.3b ; hence, photometric quality and completeness should 
be excellent for our sample. 

Stellar sources with angular separations > 0.5-0.7" are re- 
solved in SDSS photometry; we determined this empirically 
from a search in the NEIGHBORS table. For sources brighter 
than r ~ 20, the astrometric accuracy is 45 mas rms per coor- 
dinate while the re lative astrometry between filters is accurate 
to 25-35 mas rms jPier et al.ll2003l) . 

The SDSS DR7 photometric database has more than 180 
million stellar sources ( Abazaiian et al. 2009); to select a sam- 
ple o f low-mass stars, we followed the procedures outlined in 
IBIOI and required r — i > 0.3 and i — z > 0.2, which repre- 
sents the locus of sample K5 or later dwarfs. The STAR table 
was used to ensure that all of the selected objects had mor- 
phology consistent with being point sources (TYPE = 6) and 



were not duplicate detections of the same source (PRIMARY 
This yielded a sample of > 109 million low-mass stars with 
r ~ 14-24, at distances of ^0.01-5 kpc from the Sun. Fig- 
ure[T]shows a graphical flow chart of the selection process and 
criteria used; the steps needed in identifying kinematic com- 
panions, which are detailed in this paper, are shown in gray 
boxes. The white boxes show subsets where kinematic infor- 
mation is not available; however, it is still possible to identify 
binary companions based on their close proximity. We will 
discuss the selection of binaries, without available kinematic 
information, in a future paper. 

2.2. Proper Motions 

The SDSS/USNO-B matched catalog jMunn et al J 120041) . 
which is integrated into the SDSS database in the Proper- 
Motions table, was used to obtain proper motions for this 
study. We used the proper motions from the DR7 catalog; 
the earlier data releases had a systematic error i n the calcu- 
lation of proper motion in right ascension (see iMunn et al] 
I2OOI for details). This catalog uses SDSS galaxies to recali- 
brate the USNO-B positions and USNO-B stellar astrometry 
as an additional epoch for improved proper motion measure- 
ments. The ProperMotions catalog is resolution Hmited 
in the USNO-B observations to 7" and is 90% complete to 
g = 19. "fl corresponding to the faintness limit of the POSS- 
II plates used in USNO-B. The completeness also drops with 
increasing proper motion; for the range of p roper motions in 
our sample (^ =40-350 mas yr~^; see i) |2.3| below), the com- 
pleteness is ~ 85%. The typical 1 a error is 2.5-5 mas yr~^ 
for each component 

2.3. Quality Cuts 

To ensure that the resultant sample of binaries is not con- 
taminated due to bad or suspect data, we made a series of cuts 
on the stellar photometry and proper motions. With > 109 
million low-mass stars, we could afford to be very conserva- 
tive in our quality cuts and still have a reasonable number of 
stars in our input sample. We restricted our sample, as shown 
in Figure [T] to stars brighter than r = 20 and made a cut on 
the standard quality flags on the riz magnitudes (PEAKCEN- 
ter, notchecked, psf_flux_interp, interp_center, 
BAD_COUNTS_ERROR, SATURATED - all of which are re- 
quired to be Ofl which are the only bands pertinent to low- 
mass stars and the only b ands used in our an alysis. 

On the proper motions. iMunn et aP (120041 ) recommended a 
minimum total proper motion of 20 mas yr~^ and cuts based 
on different flags for a "clean" and reliable sample of stellar 
sources. Therefore, we required that each star (i) matched an 
unique USNO-B source within 1" (MATCH > 0), (ii) had no 
other SDSS source brighter than g = 22 within 7" (DIST22 
> 7), (iii) was detected on at least 4 of the 5 USNO-B plates 
and in SDSS (nFit = 6 or (nFit = 5 and (O < 2 or J < 2))), 
(iv) had a good least-squares fit to its proper motion (sigRA 
< 1000 and SIGDec < 1000), and (v) had 1 a error for both 
components less than 10 mas yr^^. 

A challenge inherent in using a deep survey like SDSS 
to identify CPM binaries, is that most of the stars are far 
away and, therefore, have small proper motions. To avoid 

^ We note that binaries separated by less than '^1" might appear elongated 
in SDSS photometry and might not be listed in the STAR table. 

* This coiTesponds to r = 18.75 for K5 and r = 18.13 for M6 dwarfs. 

^ A PRIMARY object is already selected to not be BRIGHT and NODE- 
BLEND or DEBLEND_NOPEAK 
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Figure 1. A graphical representation of the selection process used in the identification of SLoWPoKES pairs. The gray boxes show the steps involved in 
SLoWPoKES, with the results presented in this paper. The pairs without proper motions (shown in white boxes) will be presented in a follow-up paper. 



confusing real binaries with chance alignments of stars at 
large distances, where proper motions are similar but small, 
a minimum proper motion cut needs to be applied. Figure |2] 
shows the distribution of candidate binaries (selected as de- 
tailed in § 13.1b with minimum proper motion cuts of 20, 30, 
40, and 50 mas yr^^; the histograms have been normalized 
by the area of the histogram with the largest area to allow for 
relative comparisons. All four distributions have a peak at 
small separations of mostly real binaries, but the proportion 
of chance alignments at wider separations becomes larger and 
more dominant at smaller /i cutoffs. 

If our aim were to identify a complete a sample of binaries, 
we might have chosen a low fj, cutoff and accepted a relatively 
high contamination of false pairs. S uch wa s the approach in 
the recent study of SDSS binaries by lSIJOSi who used fi > 15 
mas yr^^. Our aim, in this paper, was to produce a "pure" 
sample with a high yield of bona fide binaries. Thus, in our 
search for CPM pairs we adopted a minimum proper motion 



of /i = 40 mas yr^^ for our low-mass stellar sample since the 
number of matched pairs clearly declines with increasing A9 
with this cutoff (Figure|2]i. While this still allows for a number 
of (most likely) chance alignments, they do not dominate the 
sample and can be sifted more effectively as discussed below. 

Thus, on the ^15.7 million low-mass stars that satisfied the 
quality cuts, we further imposed a /i > 40 mas yr^^ cut on 
the total proper motion that limits the input sample to 577459 
stars with excellent photometry and proper motions. As Fig- 
ure [T] shows, this input sample constitutes all the stars around 
which we searched for companions. Figure [3] shows the dis- 
tributions of photometric distances, total proper motions, r — i 
vs. i — z color-color diagram, and the r vs. r — z Hess diagram 
for the input sample. As seen from their r — z colors in the 
Hess diagram, the sample consists of K5-L0 dwarfs. Metal- 
poor halo subdwarfs are also clearly segregated from the disk 
dwarf population; however, due to a combination of the mag- 
nitude and color limits that were used, they are also mostly K 



SLoWPoKES 



5 



10- 



10- 



10- 




Lih rJlr I j. 



10 20 30 40 50 60 

AG (") 

Figure 2. The angulai' separation distribution of candidate pairs with mini- 
mum proper motions cut-offs of 20 (dotted), 30 (dashed), 40 (soUd), and 50 
(long dashes) mas yr^ ^ . These candidates were selected by matching an gular 
separations, proper motions, and photometric distances as described in ^ 13.11 
All histograms have been normalized by the area of the histogram with the 
largest area to allow for relative comparisons. It can be clearly deduced that a 
low proper motion cut-off will cause a large proportion of chance alignments. 
We adopted a cut-off of fJ- > 40 mas yr^ ^ as it allows for identification of 
CPM pairs with a reasonable number of false positives that are later sifted 
with our Galactic model. 

subdwarfs and are limited to only the earliest spectral types 
we probe. 

2.4. Derived properties 
2.4. 1 . Photometric Distances 

Disk dwarfs (DD): We determined the distances to the DDs in 
our sample by using photometric parallax relations, measured 
empirically with the SDSS stars. For M and L dwarfs (~M0- 
LO; 0.94 <r-z < 4.34), we used the relation derived bv lBld 
based on data from D. Golimowski et al. (2010, in prep.). For 
higher mass dwarfs (--05-K9; -0.72 <r - z < 0.94), we fit 
a third-order polynomial to the data reported in ICovey et"aD 
dlOOTl) : 

Mr = 3.14+7.54 (r-z)-3.23(r-z)2+0. 58 (r-z)3. (1) 

In both of the above relations, we used extinction-corrected 
magnitudes. Ideally, Mr would have been a function of 
both color and metallicity; but the effect of metallicity is not 
quantitatively known for low-mass stars. This effect, along 
with unresolved binaries and the intrinsic width of the MS, 
cause a non-Gaussian scatter of ^ 0.3-0.4 magnitudes in the 
photometric parall ax relations (IWest. Walkowicz. & Hawleyl 
l200l[s!j08l: lBToh. Since we are matching the photometric 
distances, using the smaller error bars ensures fewer false 
matches. Hence, we adopted 0.3 magnitudes as our error, 
implying a 1 tr error of ~14% in the calculated photometric 
distances. 

Subdwarfs (SD): Reliable photometric distance relations are 
not available for SDs, so instead we used the relations for 
DDs above. As a result of appearing under-luminous at a 
given color, their absolute distances will be overestimated. 
However, the relative distance between two stars in a physical 
binary should have a small uncertainty. Because we are 
interested in determining if two candidate stars occupy the 
same volume in space, the relative distances will suffice. 
While photometric parallax relations based on a few stars 
of a range of metallicities (iReid.1998; .Reid et al...2001bi) or 



caUbrated for solar-type stars ( llvezic et al.l l2008h do exist 
and can provide approximate distances for low-mass SDs, 
we refrained from using them due to the large uncertainties 
involved. 

White dwarfs (WD): We calculated the photometric 
distances to WDs using the algorithm used by [Harris et"an 
(120061) : ugriz magnitudes and the u — g, g — r, r — i, and i — z 
colors, cor rected for extinction, were fitted to th e WD c ooling 
models of [Bergeron. Wesemael. & Beaucharnpl ( Il995h to get 
the bolometric luminosities and, hence, the distance^. The 
bolometric luminosity of WDs are a function of gravity as 
well as the composition of its atmosphere (hydrogen/helium), 
neither of which can be determined from the photometry. 
Therefore, we assumed pure hydrogen atmospheres with 
log g = 8.0 to calculate the distances to the WDs. As a 
result, distances derived for WDs with unusually low mass 
and gravity (^10% of all WDs) or with unusually high mass 
and gravity ('^15% of all WDs) will have larger uncertainties. 
Helium WDs redder than g — i > 0.3 will also have discrepant 
distances ([Harris et al.ll2006h . 

2.4.2. Spectral Type & Mass 

The spectral type of all (05-M9; -0.6 <r - z <4.5) disk 
dwarfs and subdwarfs were inferred from their r — z colors 
using the following two-part fourth-order polynomial: 

^2 , 



SpT = 



50.4 + 9.04 (r - z) - 2.97 (r - 

0.516 (r - zf - 0.028 (r - z)*; for 05-K5^ 
33.6 + 47.8 (r - z) - 20.1 (r - z f- 

33.7 (r - zf - 58.2 (r - z)^; for K5-M9 



(2) 



where SpT ranges from 0-67 (00=0 and M9=67 with all 
spectral types, except K for which K6, K8, and K9 are not de- 
fined, having 10 su btyp es) and is based on the data reported in 
[Covey et al.r([2007[) and [West etal] ([20081) . The spectral types 
should be correct to ±1 subtype. 

Similarly, the mass of B8-M9 disk dwarfs and subdwarfs 
were determined from their r — z colors using a two-part 
fourth-order polynomial,: 



1.21 - 1.12 (r - z) + 1.06 (r - z)^- 



2.97 (r 



z)'i;forB8-K5 
0.362 (r-z)- 0.539 (r-z)2+ ''^> 
zf - 0.016 (r - z)4; for K5-M9 



2.73 (r 
0.640 

0.170 (r - zf - 0.016 (r 

based on the data reported in [Kraus & Hillenbrand! ( [2007h . 
who used theoretical models, supplemented with observa- 
tional constrains when needed, to get mass as a function of 
spectral type. The scatter of the fit, as defined by the median 
absolute deviation, is ^2%. 

As the polynomials for both the spectral type and mass are 
monotonic functions of r — z over the entire range, the com- 
ponent with the bluer r — z color was classified as the primary 
star of each binary found. 

3. METHOD 

3.1. Binary Selection 

Components of a gravitationally bound system are expected 
to occupy the same spatial volume, described by their semi- 
axes, and to move with a common space velocity. To identify 

^ SDSS ugriz magnitudes and colors for the WD 
cooling models are available on P. Bergeron's website: 
|http : / /www .astro ■ umontreal ■ ca/ -bergeron/CoolingModels7] 
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Figure 3. The characteristics of the sample of 577459 low-mass stars in the SDSS DR7 photometric catalog that forms our input sample (clockwise from top 
left): the photometric distances, the total proper motions, r vs. r — 2 Hess diagram, and the r — i vs. i — z color-color plot. In addition to quality cuts on 
proper motion and photometry, we require that stars in this sample to be low-mass (r — i > 0.3 and i — z > 0.2) and have a relatively high proper motion (p > 
40 mas yr~^). The contours show densities of 10-3000 in the color-color plot and 10-300 in the Hess diagram. 



physical binaries in our low-mass sample, we implemented a 
statistical matching of positional astrometry (right ascension, 
a, and declination, 6), proper motion components and 
l-is), and photometric distances (d). The matching of distances 
is an improve ment to the methods of previous searches for 
CPM doubles jLepine & BongiornoliOOTtlChaname & Gouldl 
12001 lHalbwachslll986h and serves to provide further confi- 
dence in the binarity of identified systems. 

The angular separation, A9, between two nearby point 
sources, A and B, on the sky can be calculated using the small 
angle approximation: 



AO ~ y/iaA - as) 2 cos (5a COS (5b + {5a - Sb) 



(4) 



We searched around each star in the input sample for all stel- 
lar neighbors, brighter than r = 20 and with good photom- 
etry and proper motions, within A9 < 180" in the SDSS 
photometric database using the Catalog Archive Server query 
too0 Although CPM binaries have been found at much larger 
angular s eparations (up to 900" in [ Chaname" .fe Gouldl I2OOI 
1500" in iLepine & Bongiornd I^OOtT 570" in lFahertveTan 
I2OIOI) . the contamination rate of chance alignments at such 
large angular separations was unacceptably high in the deep 
SDSS imaging (see Figure|2]i. In addition, searching the large 
number of matches at larger separations required large com- 
putational resources. However, since the SDSS low-mass star 
sample spans considerably larger distances than in previous 
studies (see below), our cutoff of 180" angular separation 
probes similar physical separations of up to '^0.5 pc, which 
is c omparable to the typica l size of prestellar cores (0.35 
pc: iBenson & MversI [l989t IClemens. Yun. & Heyeil [l99Tt 
iJessop & Ward-Thompsonl2000^ . 

^ |http : //cas jobs ■ sdss ■ org/Cas Jobs/| 



For all pairs that were found with angular separations of 
7" < A9 < 180", we required the photometric distances to 
be within 

Ad < (1 aAd, 100 pc) (5) 
and the proper motions to be within 



AHa 



< 2 



(6) 



where Ad, AyUc and Ai^ls are the scalar differences between 
the two components with their uncertainties calculated by 
adding the individual uncertainties in quadrature. An absolute 
upper limit on Ad of 100 pc was imposed to avoid Aii being 
arbitrarily large at the very large distances probed by SDSS; 
hence, at d > 720 pc, distances are matched to be within 
100 pc and results in a significantly lower number of candi- 
date pairs. The proper motions are matched in 2-dimensional 
vector space, instead of just matching the total (scalar) proper 
motion as has been frequently done in the past. The latter 
approach allows for a significant number of false positives as 
stars with proper motions with the same magnitude but differ- 
ent directions can be misidentified as CPM pairs. 

For our sample, the uncertainties in proper motions are 
almost always larger than the largest possible Keplerian or- 
bital motions of the identified pairs. For example, a binary 
with a separation of 5000 AU and Mtot= 1 M© at 200 
pc, which is a typical pair in the resultant SLoWPoKES sam- 
ple, will has a maximum orbital motion of ~0.27 mas yr^^, 
much smaller than error in our proper motions (typically 2.5- 
5 mas yr^^). However, a pair with a separation of 500 AU 
pair and Mtot= 1 Mq at 50 pc has a maximum orbital mo- 
tion of ~^5.62 mas yr~^, comparable to the largest errors in 
component proper motions. Hence, our algorithm will reject 
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Figure 4. The spatial distribution of SLoWPoKES binaries, sliown in equa- 
torial coordinates, closely follows the SDSS footprint. The upper and right 
panels show the histogram of Right Ascension and DecHnation in 5° bins. 
We rejected the 118 CPM candidates found in the direction of open clus- 
ter IC 5146 (shown in red pluses and dashed histograms) due to the highly 
anomalous concent ratio n of pairs and possibly contaminated photometry due 
to nebulosity (see ii l3.U . 

such nearby, relatively tight binaries. 

From the mass estimates and the angular separations from 
the resulting sample, we calculated the maximum Keplerian 
orbital velocities, which are typically less than 1 mas yr^^; 
only 104, 34, and 3 out of a total of 1342 systems exceeded 
1, 2, and 5 mas yr^^, respectively. More importantly, only 
7 pairs had maximum orbital velocities greater than l-cr error 
in the proper motion; so apart from the nearest and/or tightest 
pairs, our search should not have been affected by our restric- 
tions on the proper motion matching in Eq. (|6]l. 

Applying these selection criteria to the stellar sample de- 
scribed in § 12.11 we found a total of 3774 wide CPM binary 
candidates, where each has at least one low-mass component. 
Among these, we found 118 pairs shown in red pluses and 
dashed histograms in Figure |4l concentrated in a 20' x 150' 
stripe, in the direction of open cluster IC 5 146. As there is sig- 
nificant nebulosity that is not reflected in the extinction values, 
the photometry and, hence, the calculated photometric dis- 
tances are not reliable. In addition, while the kinematics are 
not characteristic of IC 5146, they are more likely to be part 
of a moving group rather than individual CPM pairs. Thus, 
we rejected all of these candidates. No other distinct struc- 
tures, in space and kinematics, were found. Then, we made 
the quality cuts described in §[23] on the companions; 906 
and 1085 companions did not meet our threshold for the pho- 
tometry and proper motions, respectively, and were rejected. 
The majority of these rejections were near the Galactic Plane, 
which was expected due to higher stellar density. Thus, at the 
end, the resulting sample had 1598 CPM double candidates 
from the statistical matching (Figure[T]i. Inherent in statistical 
samples are false positives, arising from chance alignments 
within the uncertainties of the selection criteria. For any star, 
the probability of chance alignment grows with the separation 
making the wider companions much more likely to be chance 
alignments despite the selection criteria we just implemented. 
Hence, it is necessary to complete a detailed analysis of the 
fidelity of the sample. 

3.2. Galactic Model: Assessing False Positives in the Binary 

Sample 



To quantitatively assess the fidelity of each binary in our 
sample, we created a Monte Carlo based Galactic model that 
mimics the spatial and kinematic stellar distributions of the 
Milky Way and calculates the likelihood that a given binary 
could arise by chance from a random alignment of stars. 
Clearly, for this to work the underlying distribution needs to 
be carefully constructed such that, as an ensemble, it reflects 
the statistical properties of the true distribution. The model 
needs to take into account the changes in the Galactic stellar 
density and space velocity over the large range of galacto- 
centric distances and heights above/below the Galactic disk 
plane probed by our sample. Previous studies, focused on 
nearby binaries, have been able to treat the underlying Galac- 
tic stellar distribution as a simpl e random distribution in two- 
dimensional space. For example. iLepine & Bongiornol(l2007h 
assigned a random shift of 1-5° in Right Ascension for the 
secondary of their candidate pairs and compared the resultant 
distribution with the real one. As they noted, the shift cannot 
be arbitrarily large and needs to be within regions of the sky 
with similar densities and proper motion systematics. With 
stars at much larger heliocentric distances in our sample, even 
a 1° shift would correspond to a large shift in Galactic position 
(1° in the sky at 1000 pc represents ^17.5 pc). Therefore, we 
could not use a similar approach to assess the false positives 
in our catalog. 

The Galactic model is built on the canonical view that the 
Galaxy comprises three distinct components — the thin disk, 
the thick disk, and the halo — that can be cleanly segregated 
by their age, metallicity, a nd kinematics (Bahcall & Soneira] 
ll98d:lGilmore. Wyse. & Kuijken 1 989; Majewski 1993) . We 
note that some recent work has argued for the disk to be a 
continuum instead of two distinct components jlvezic et al.l 
120081) or for the halo to be compo sed of two distinct com- 
ponents (ICaroUo et al.ll2008l l2009l) . However, for our pur- 
poses, the canonical three-component model was sufficient. 
We also did not try to model the over-densities or under- 
densities, in positional or kinematic space, caused by co- 
moving groups, open clusters, star-forming regions, or Galac- 
tic streams. If such substructures were found in t he S DSS 
data, they were removed from our sample (see § 13.1b . In 
essence, this model strictly describes stars in the field and 
produces the three-dimensional position and two-dimensional 
proper motion, analogous to what is available for the SDSS 
photometric catalog. 

3.2. 1 . Galactic stellar density profile 

In the canonical Galactic model, the stellar densities (p) of 
the thin and the thick disks, in standard Galactic coordinates 
R (Galactic radius) and Z (Galactic height), are given by 

,(i?,Z) = p (i?0,O) e «th.„ e ^tM„ 



Pthi: 



(7) 
(8) 



where H and L represent the scale height above (and below) 
the plane and the scale length within the plane, respectively. 
The halo is described by a bi-axial power-law ellipsoid 



Phalo(i?, •Z') = p (i?0,O) 



(9) 



with a halo flattening parameter q and a halo density gradient 
Thaio- The three profiles are added together, with the appro- 
priate scaUng factors, /, to give the stellar density profile of 
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Table 1 

Galactic Structure Parameters 



Table 2 

Galactic Kinematics 



Component 


Parameter name 


Parameter description 


Adopted Value 




p{Re,0) 


stellar density 


0.0064 




/thin 


fraction^ 


"J thick "J halo 


thin disk 


-Wthin 


scale height 


260 pc 






scale length 


2500 pc 




/thick 


fraction" 


9% 


thick disk 


-Wthick 


scale height 


900 pc 




-^'thick 


scale length 


3500 pc 




/halo 


fraction" 


0.25% 


halo 


''halo 


density gradient 


2.77 




g (= c/a)" 


flattening parameter 


0.64 


Note. 


The parameters 


were measured using M 


dwarfs for the 



disk (Bochanski et al. 2010) and main-sequence turn-off stars for the halo 
Ouric et aU200a) in the SDSS footprint. 
" Evaluated in the solar neighborhood 



' Assuming a bi-axial ellipsoid with axes a and c 



the Galaxy: 

P{R,Z): 



/thin Pthin + /thick Pthick + /halo Phalo- (10) 



The scaling factors are normalized such that /thin + /thick + 
/halo = 1- With the large number of stars imaged in the 
SDSS, robust stellar density functions have been m easured for 
the thin an d thick disks using the low-mass stars ( lJuric et al.l 
120081: IbToI) and for the halo using he main-sequence turn-off 
stars (lJuric et al.ll2008h . The values measured for the disk in 
the two studies are in rough agreement. We adopt ed the disk 
parameters fromlBlQl and the halo parameters from lJuric et alJ 
d2008li) : Table[T]summarizes the adopted values. 



3.2.2. Galactic Kinematics 

Compared to the positions, the kinematics of the stellar 
components are not as well characterized; in fact, apart from 
their large velocity dispersions, little is known about the halo 
kinematics. We seek to compare the proper motions of a 
candidate pair with the expected proper motions for that pair 
given its Galactic position. Thus, we found it prudent to (i) ig- 
nore the halo component, with its unconstrained kinematics, 
at distances where its contributions are expected to be mini- 
mal and (ii) limit our model to a distance of 2500 pc, which 
corresponds to the Galactic height wh ere the number of halo 
stars begins to outnumber disk stars ( lJuric e t al. 20QS). In 
practice, all the SLoWPoKES CPM pairs, with the exceptions 
of subdwarfs for which distances were known to be overesti- 
mated, were within ~1200 pc (see Figure |9]below); so we did 
not introduce any significant systematics with these restric- 
tions. 

An ensemble of stars in the Galactic Plane can be character- 
ized as having a purely circular motion with a velocity, Vdic- 
The orbits become more elliptic and eccentric over time due 
to kinematic heating causing the azimuthal velocity, V^, to 
decrease with the Galactic height, Z. However, the mean ra- 
dial (Vr) and perpendicular (Vz) velocities for the ensemble 
at any Z remains zero, with a given dispersion, as there is no 
net flow of stars in either direction. In addition, this random- 
ization of orbits also causes the asymmetric drift, Va, which 
increases with the age of stellar population and is equivalent 
to '^lO km s^^ for M dwarfs. Hence, the velocities of stars in 
the Galactic disk can be summarized, in Galactic cylindrical 



Galactic component Velocity 


k 


n 




U 


7.09 


0.28 


thin disk 


V 


3.20 


0.35 




W 


3.70 


0.31 




U 


10.38 


0.29 


thick disk 


V 


1.11 


0.63 




W 


0.31 


0.31 



Note. — The constants in the power law, 
<y(Z) = k\Z\", that describes the velocity 
dispersions of the stars in the thin and thick 
disks. The velocity dispersions were measured 
from a spectroscopic s ample of low-mass dwarfs 
JBochanski et al.l2007D . 



coordinates, as: 



{V^ {Z)} = Vcn-c - - f{Z) 



(11) 



where Vdrc = 220 km s-\ Va = 10 km s-\ f{Z) = 
0.013 \Z\ - 1.56 X IQ -^ \Z\^ wa s deriv ed by fitting a poly- 
nomial to the data in IWest et all (l2008h . and Z is in par- 
sees. This formulation of V^{Z) is consistent with a stel- 
lar population composed of a faster thin disk ((V^) = 210 
km s~^) and a slower thick disk ((V^) = 180 km s^^). 
Then, we converted these galactocentric polar velocities to 
the heliocentric, Cartesian UVW velocities. The UVW veloc- 
ities, when complemented with the dispersions, can be con- 
ver ted to a two-dimensional prop er motion (and radial veloc- 
ity; |johnson^^odeAlo3[T987b, analogous to our input cat- 
alog. We used the UV W velocity dispersions measured for 
SDSS low-mass dwarfs (iBochanski et al.ll2007h . All the dis- 
persions were well described by the power-law 



a{Z) = k \Z\ 



(12) 



where the values of constants k and n are summarized in Ta- 
ble 12] As the velocity dispersions in IBochanski et al.l (120071) 
extend only up to ^1200 pc, we extrapolated the above equa- 
tion for larger distances. While the velocity ellipsoids of F 
and G dwarfs ha ve been measured to larger distances (e.g. 
iBond et akl 120091) . we preferred to use the values measured 
for M dwarfs for our low-mass sample. 



3.2.3. The Model 

By definition, a chance alignment occurs because two phys- 
ically unassociated stars randomly happen to be close together 
in our line of sight (LOS), within the errors of our measure- 
ments. Due to the random nature of these chance alignments, 
it is not sufficient to estimate the probability of chance align- 
ment along a given LOS simply by integrating Eq. (ITOl i. This 
would tend to underestimate the true number of chance align- 
ments because the density profiles in Eqs. (jTHUl are smooth 
functions that, in themselves, do not include the random scat- 
ter about the mean relation that real stars in the real Galaxy 
have. Thus, the stars need to be randomly redistributed spa- 
tially about the average given by Eq. (ITOl i in order to properly 
account for small, random fluctuations in position and veloc- 
ity that could give rise to false binaries in our data. 
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In principle, one could simulate the whole Galaxy in this 
fashion in order to determine the probabilities of chance align- 
ments as a function of LOS. In practice, this requires exorbi- 
tant amounts of computational time and memory. Since our 
aim was to calculate the likelihood of a false positive along 
specific LOSs, we, instead, generated stars in much smaller 
regions of space, corresponding to the specific binaries in our 
sample. For example, a 30' x 30' cone integrated out to a 
distance of 2500 pc from the Sun will contain at least a few 
thousand stars within any of the specific LOSs in our sample. 
The number of stars is large enough to allow for density vari- 
ations similar to that of the MiUcy Way while small enough to 
be simulated with ease. With a sufficient number of Monte 
Carlo realizations, the random density fluctuations along each 
LOS can be simulated. We found that 10^ realizations allowed 
for the results to converge, within ^-^0.5%. 

We implemented the following recipe to assess the false- 
positive Ukelihood for each candidate pair using our Galactic 
model: 

(i) The total number of stars in the LOS volume defined 
by a 30' x 30' area, centered on the a and (5 of a given 
binary, over heliocentric distances of 0-2500 pc was 
calculated by integrating Eq. ( fTOl i in 5 pc deep, discrete 
cylindrical "cells." 

Integrating Eq. ^ for A9 = 30' and d = 0-2500 pc 
resulted in ^3300-1580 stars, with the higher numbers 
more likely to be the LOSs along the Galactic Plane. 
This number of stars, when randomly redistributed in 
the entire volume, was more than enough to recreate 
over-densities and under-densities. 

(ii) The stars were then distributed in three-dimensional 
space defined by the LOS using the rejection method 
(iPress et al.ll 19921) . generating a, 6, and d for each star 
The rejection method ensured that the stars were ran- 
domly distributed while following the overlying distri- 
bution function, which, in this case, was the stellar den- 
sity profile given by Eq. (fTOl i. 

The model did an excellent job of replicating the ac- 
tual distribution of stars in the three-dimensional space. 
The red histograms in Figure |5] show the number of 
stars from the center of the LOS as a function of an- 
gular separation for the model (dashed lines) and the 
data (solid lines), averaged over a ll LO Ss where candi- 
date binaries were identified in § 13.11 There is an ex- 
cess of pairs at close separations, a signature of gen- 
uine, physically associated pairs, while the two distri- 
butions follow the same functional form at larger sep- 
arations, where chance alignments dominate. The in- 
creasing number of pairs with angular separation is an 
evidence of larger volume that is being searched. In- 
tegrating the model predicts ^8.0 stars within a search 
radius of 60" and ^71 within 180". In other words, for 
a typical LOS with candidate binaries, we would, on 
average, expect to find 8 chance alignments within 60" 
and ^71 within 180" when considering only angular 
separation. 

For the 1 80" radius around each LOS, the average num- 
ber of stars in the model and the data were within a few 
a of each other, with the largest deviations found along 
the LOSs at low Galactic latitudes or large distances. 
As our model only integrated to d — 2500 pc, devia- 
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Figure 5. The distribution of the (average) number of stars found around the 
LOSs around our binary candidates as calculated from the SDSS DRV data 
(solid histograms) and our Galactic model (dashed histograms). All optical 
pairs (red), pairs with matching distances (blue), and pairs with matching 
distances and proper motion components (puiple) are shown. The excess at 
small separations, a signature of real pairs, is enhanced as additional proper- 
ties are matched. Unlike the rest of the paper, we counted all stars with r < 
22.5, without any quality cuts, for this plot to do a realistic comparison with 
the model. Note that (A^) denotes the number of stars found at that angular 
separation (counted in 1" bins). We used the raw parameters for the Galactic 
model tJuric et al. 2008. BlQ). As a whole, the model does an excellent job 
of mimicking the number of stars and their spatial and kinematic distribution 
in the Milky Way along typical SDSS LOSs. 

tions for LOSs at large distances were expected. Devi- 
ations for LOSs at low Galactic latitudes are reasonable 
as the parameters in Eqs. ^ and ^ as not as well con- 
strained along the Galactic disk. Hence, we concluded 
that the rejection method used in redistributing the stars 
in the 30' x 30' LOS is correctly implemented and that 
the model successfully replicates the three dimensional 
spatial distribution of the stars in the Galaxy. 

When the distances were matched, the number of pairs 
decreased as chance alignments were rejected (blue his- 
tograms in Figure|5]i; there are, on average, only ^0.41 
and ^3.6 chance alignments within 60" and 180" , re- 
spectively. Note that matching distances, in addition to 
the angular separation, significantly enhanced the peak 
at the small separations. 

(iii) Based on the Galactic position of each randomly gener- 
ated star, mean UVW space velocities and their disper- 
sions were generated based on Eq. ( fTTT i and Eq. ( fT2b , 
respectively. Proper motions were then calculated 
by using the inverse of th e algorithm outlined in 
iJohnson & SoderblomI (11987'). These generated proper 
motions represent the expected kinematics of stars at 
the given Galactic position. 

Figure|6]shows the comparison between the proper mo- 
tions in the SDSS/USNO-B catalog and our model; 
component proper motions of stars within 60' of all 
LOSs, where candidate pairs were found. For the pur- 
poses of this plot, we restricted the stars to be within 
1200 pc and of spectral type K5 or later so we could 
compare kinematics with a sample similar to our resul- 
tant catalog. We also did not compare the distributions 
at ji a or: jig < 10 mas yr~^ where the proper motions 
are comparable to the l-cr errors and, hence, not reli- 
able. Since our initial sample was has /i > 40 mas yr^^, 
stars with the small proper motion component are either 
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Figure 6. Proper motion distributions for stars within 60' of LOSs along all 
identified candidate pairs in SDSS/USNO-B (solid histograms) and in our 
Galactic model (dashed histograms). For the purposes of this plot, we re- 
stricted the stars to be of spectral types K5 or later and be within 1200 pc to 
avoid being skewed by systematic differences. In addition, we do not com- 
pare proper motion components < 10 mas yr^^ as they comparable to the 
1-(T errors and, hence, not reliable (see text). Again, the kinematics of the 
thin and thick disk of the Milky Way are very well reproduced by our Monte 
Carlo model. 

rejected or their motion is dominated by the other com- 
ponent in our analysis. As evidenced in the figure, our 
model reproduced the overall kinematic structure of the 
thin and thick disks. 

When the proper motions were matched (purple his- 
tograms in Figure |5]l, in addition to the angular sepa- 
ration and distances, the number of chance alignments 
fell drastically, especially at the smaller separations. In 
fact, at /S.9 < 15", there were < 10"'* chance align- 
ments, on average; even at A0= 180", the real pairs 
outnumbered the chance alignments by a factor of 4- 
5. Cumulatively, for a typical LOS, there were, on 
average, only ^0.0097 and ^0.086 chance alignments 
within 60" and 180", respectively. As a result of match- 
ing distance and proper motions components, the num- 
ber of chance alignments within 180" were reduced by 
a factor of --800. 

(iv) In the model galaxy, we repeated the selection process 
used to find CPM pair s in the SDSS photometric cat- 
alog, as described in § 13.11 To avoid double-counting, 
the input coordinates of the LOS were considered to be 
the primary star Note that we did not intend to model 
and reproduce both stars of a given pair but wanted to 
see if a random chance alignment could produce a com- 
panion for a given primary. Hence, each additional star 
that was found to satisfy Eqs. (|5]l and ^ was counted 
as an optical companion. The average number of com- 
panions found in the 10^ Monte Carlo realizations is the 
probability of chance alignment or the probability that 
the candidate pair is a false positive, Pf, for that candi- 
date pair The number of realizations sets the resolution 
ofPf at IQ-^. 

(v) Finally, this was rep eated for all candidate pairs that 
were found in § 13.11 Figure [7] shows the distribution 
of Pf all of the candidate pairs; we further discuss this 
distribution is i) l3.3l 

To conclude, the result of our Galactic model was a 30' x 30' 
area of the sky, centered around the given CPM pair, with the 
surrounding stars following the Galactic spatial and kinematic 
distributions of the Milky Way. Each star in this model galaxy 
was described by its position (a, 5, and d) and proper motion 



(/XQ.and /i^), same as is available for SDSS photometric cata- 
log. Based on the above results, we concluded that the Galac- 
tic model sufficiently reproduced the five-dimensional (three 
spatial and two kinematicFl) distribution of stars along typical 
LOSs in the Milky Way and, thus, allowed for the calculation 
of probability that a given CPM double is a chance alignment. 

3.3. Fidelity 

As described in the previous section, we have implemented 
a very stringent selection algorithm in identifying the CPM 
pairs. In addition to only including objects with the most ro- 
bust photometry and proper motions, we also used a relatively 
high proper motion cut of /i > 40 mas yr^* for the input sam- 
ple, which considerably decreased the number of low-mass 
stars. As described above, an algorithm optimized to reject 
the most false positives, even at the expense of real pairs, was 
used in the statistical matching of angular separation, photo- 
metric distance, and proper motion components. Lastly, we 
used the Galactic model to quantify the probability of chance 
alignment, Pf , for each of the 1558 candidate pairs. 

Normally all candidate pairs with Pf < 0.5, i.e. a higher 
chance of being a real rather than a fake pair, could be used 
to identify the binaries. However, to minimize the number of 
spurious pairs, we required 

Ff < 0.05. (13) 

for a pair to be classified as real. Here we note that Pf rep- 
resents the false-alarm probability that a candidate pair iden- 
tified by matching angular separation, photometric distance, 
and proper motion components, as described by Eqs. (|5]l and 
is a real pair; it is not the probability that a random low- 
mass star is part of a wide binary. 

Making the above cut on Pf resulted in a catalog of 1342 
pairs, with a maximum of 5% or 67 of the pairs expected to 
be false positives. However, as a large number of the pairs 
have extremely small Pf (see Figure |2ll, the number of false 
positives is likely to be much smaller Adding up the Pf for 
the pairs included in the catalog gives an estimated 22 (1.65%) 
false positives. In other words, the overall fidelity of SLoW- 
PoKES is 98.35%. This is a remarkably small proportion for 
a sample of very wide pairs, especially since they span a large 
range in heliocentric distances and is a testament to our selec- 
tion criteria. For example, if the proper motion components 
were matched to within 2 cr, Eq. dljl l would have rejected > 
60% of the candidates. Our choice of Eq. (fT3] i is a matter of 
preference; if a more efficient (or larger) sample is needed, it 
can be changed to suit the purpose. 

To get a first-order approximation of how many real bi- 
naries we are missing or rejecting, we applied our selec- 
tion algorithm to the rNLTT CPM catal og with 1147 pairs 
(?Cha name & Gou ld 2004). Note that IChaname & Gouldl 
(120041) did not match distances of the components, as they did 
not have reliable distance estimates available and matched to- 
tal proper motions instead of a two-dimensional vector match- 
ing in our approach. Out of the 307 rNLTT pairs, which have 
/S.9 < 180" and are within the SDSS footprint, we recover 
both components of 194 systems (63%) and one component 
of another 56 (18%) systems, within 2" of the rNLTT coordi- 
nates. Of the 194 pairs for which both components had SDSS 
counterparts, 59 (30%) pairs satisfied our criteria for proper 

The model also predicts radial velocities. We have an observational 
program underway to obtain radial velocities of the binaries for further re- 
finement of the sample. 
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Figure 7. The probability of chance alignment (Pf ) for candidate SLoW- 
PoKES pairs based on their positions and kinematics, as calculated by our 
Galactic model. We adopted Pf = 0.05 (dotted line) as our threshold for 
inclusion in the SLoWPoKES catalog. The resolution of the Monte Carlo 
simulation is 10^^, which causes the peak at log Pf = 10"''. 

motions, Eq. while 19 (10%) pairs satisfied our criteria for 
both proper motions and distances, Eqs. Q and (|6]l. In other 
words with our selection criteria, we recovered only 10% of 
the rNLTT pairs, with SDSS counterparts, as real binaries. If 
we relax our selection criteria to match within 2-cr, the num- 
ber of recovered pairs increases to 82 with matching proper 
motions and 37 with matching proper motions and distances. 
Of course, as the rNLTT is a nearby, high proper motion cat- 
alog, [Chaname & Gould (2004) had to allow for larger differ- 
ences in proper motions due to the stars' orbital motio n, w hich 
is not applicable to the SLoWPoKES sample (see §[TT]i. In 
conclusion, compared to previous catalogs of CPM doubles, 
we find a small fraction of previously identified very wide bi- 
naries. The reasons for the low recovery rate are two-fold: 

(i) the restrictive nature of our matching algorithm that rejects 
the most false positives, even at the expense of real pairs and 

(ii) improvement in the identification method — e.g. matching 
proper motions in vector space and being able to use photo- 
metric distance as an additional criterion. 



4. CHARACTERISTICS OF THE SLoWPoKES CATALOG 

Using statistical matching of angular separation, photomet- 
ric distance, and proper motion components, we have iden- 
tified 1598 very wide, CPM double candidates from SDSS 
DR7. We have built a Galactic model, based on empirical 
measurements of the stellar density profile and kinematics, to 
quantitatively evaluate the probability that each of those can- 
didate doubles are real (Figure |7]). Using Eq. (T3[ . we clas- 
sify 1342 pairs as real, associated pairs. In deference to their 
extremely slow movement around each other, we dub the re- 
sulting catalog SLoWPoKES for Sloan Low-mass Wide Pairs 
of Kinematically Equivalent Stars. Table|4]lists the properties 
of the identified pairs. The full catalog is publicly available 
on the world wide wetP^. SLoWPoKES is intended to be a 
"live" catalog of wide, low-mass pairs, i.e., it will be updated 
as more pairs are identified and as follow-up photometric and 
spectroscopic data become available. 

Figure [8] shows a collage of gri composite images, 50" on 
a side, of a selection of representative SLoWPoKES systems; 

'* |http : / / www ■ vanderbilt ■ edu/astro/ slowpokes/] 



the images are from the SDSS database. In the collage high- 
mass ratio pairs (top row; mass ratio = m2 / mi < 0.5), equal- 
mass pairs (middle row; having masses within ^5% of each 
other), white dwarf-disk dwarf pairs (bottom row, left), and 
halo subdwarf pairs (bottom row, right) are shown. Table |4] 
summarizes the different types of systems in the catalog. 



In the sections that follow, we summarize various aspects 
of the systems that constitute the SLoWPoKES catalog and 
briefly examine some of the follow-up science that can be 
pursued with SLoWPoKES. We wish to emphasize that the 
SLoWPoKES catalog is intended principally as a high-fidelity 
sample of CPM doubles that can be used for a variety of 
follow-up investigations where the reliability of each object 
in the catalog is more important than sample completeness. 
Thus, we have not attempted to fully account for all sources 
of incompleteness or bias; and we intentionally have not ap- 
plied any form of statistical "corrections" to the catalog. 

We note that some of the most important sources of incom- 
pleteness in the present catalog could be at least partially over- 
come with follow-up observations. For example, the prin- 
cipal incompleteness in SLoWPoKES arises from the lack 
of proper motions for SDSS stars that were not detected in 
USNO-B. Proper motions either do not exist or are not reliable 
for stars (i) fainter than r = 20 or (ii) within 7" of a brighter 
star The first criterion currently rules out most of the mid- 
late M dwarf companions, while the latter rejects close bina- 
ries and most hierarchical higher-order systems. However, as 
the SDSS photometric and astrometric data are available for 
these systems, their multiplicity could be verified with more 
rigorous analysis or through cross-matching with other cata- 
logs. For example, by cross-matching S DSS with 2MASS, 
an M4.5-L5 binary (Radi gan etal]l2008l) and an M6-M7 bi- 
nary (jRadigan et al. 2009) have already been identified. At 
the other end of the spectrum, SDSS saturates at r w 14, re- 
sulting in saturated or unreliable photometry for the brighter 
stars. We found that more than 1000 candidate pairs were 
rejected for this reason; with reliable follow-up photometry 
these could be added as additional genuine SLoWPoKES bi- 
naries. 

Even so, a fully volume-complete sample is likely to be im- 
possible to compile over the full ranges of spectral types and 
distances spanned by our catalog. For example, our magni- 
tude limits of 14 < r < 20 imply that while we are sensi- 
tive to K5 dwarfs at ~^250-3900 pc, we are sensitive to M5 
dwarfs at ^14-180 pc. This means (i) we are entirely insen- 
sitive to pairs with the most extreme mass ratios (i.e. a K5 
paired with an M6 or M7) and (ii) we cannot directly com- 
pare the properties of K5 and M5 spectral subtypes in iden- 
tical distance ranges. In addition, as illustrated in Figure |9l 
with our 7" < A6 < 180" search radius we are sensitive to 
companions with separations of ^700-18000 AU at 100 pc 
but -7000-180000 AU at 1000 pc. Thus, it is important with 
the current catalog that statistical determinations of ensem- 
ble system properties be performed within narrowly defined 
slices of separation and distance. We do so in the follow- 
ing subsections as appropriate; but we emphasize again that 
our intent here is primarily to characterize the SLoWPoKES 
sample and will proffer any interpretive conclusions only ten- 
tatively. 

Finally, as is evident from Table IH the current SLoW- 
PoKES catalog is clearly not well suited for study of higher- 
order multiples (i.e. triples, quadruples, etc.). Identifying 
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Table 3 

Properties of SLoWPoKES pairs 



ID" 


Position (J2000) 










Photometry'' 










As 


'!"/ 


1 


«A 


ZA 






ZB 


SLW 


(deg) 










(mag) 






J0002+29 0.515839 29.475183 


0.514769 


29.470617 


16.79 (0.02) 


15.85 (0.02) 


15.33 (0.02) 


19.35 (0.02) 17.91 (0.02) 


17.17 (0.02) 


J0004-10 1.122441 -10.324043 


1.095927 


-10.296753 


18.25 (0.02) 


17.08 (0.02) 


16.46 (0.02) 


18.66(0.02) 17.42(0.02) 


16.70 (0.02) 


J0004-05 1.223125 -5.266612 


1.249632 


-5.237299 


14.85 (0.01) 


14.31 (0.01) 


14.02 (0.02) 


19.56(0.02) 18.10(0.01) 


17.31 (0.02) 


J0005-07 1.442631 -7.569930 


1.398478 


-7.569359 


17.35 (0.02) 


16.25 (0.01) 


15.65 (0.01) 


18.78 (0.02) 17.43 (0.02) 


16.69 (0.01) 


10005-1-27 1.464802 27.325805 


1.422484 


27.300282 


18.91 (0.02) 


17.63 (0.02) 


16.94 (0.01) 


19.79(0.02) 18 


.29 (0.02) 


17.50 (0.02) 


J0006-03 1.640868 -3.928988 


1.641011 


-3.926589 


17.06(0.01) 


16.01 (0.02) 


15.44(0.01) 


17.88 (0.01) 16.62(0.02) 


15.93 (0.01) 


J0006-H08 1.670973 8.454040 


1.690136 


8.498541 


19.48 (0.02) 


18.14(0.02) 


17.45 (0.02) 


19.57 (0.02) 18, 


.15 (0.01) 


17.43 (0.02) 


J0007-10 1.917002 -10.340915 


1.924510 


-10.338869 


17.43 (0.01) 


16.50 (0.01) 


16.02 (0.03) 


18.33 (0.01) 17.24(0.01) 


16.66 (0.03) 


J0008-07 2.135590 -7.992694 


2.133660 


-7.995245 


17.33 (0.01) 


16.35 (0.02) 


15.86 (0.02) 


18.14(0.01) 16 


.98 (0.02) 


16.36 (0.02) 


J0009-H15 2.268841 15.069630 


2.272453 


15.066457 


18.52 (0.02) 


17.00 (0.02) 


16.09(0.01) 


18.39(0.02) 16.83 (0.02) 


15.95 (0.01) 


Table 3 




Proper Motion 




Distance"^ 


Spectral Type'' 




Binary Information 












dA dB 


A 


B 


AO 




Ad BH 




Class' 




( mas yr^ ^ ) 




(pc) 






(") 


( mas yr^ ' ) 


(pc) (10*0 g;.gj,) 


(%) 




198 (2) 


38 (2) 197 (3) 


35 (3) 


341 301 


M1.7 


M3.6 


16.8 


2.9 


39 58.08 


0.000 


SD 


43 (3) 


-4 (3) 36 (4) 


-8 (4) 


362 338 


M2.7 


M3.1 


135.9 


7.5 


23 3.94 


0.036 


DD 


101 (2) 


11(2) 99(4) 


8(4) 


301 333 


K7.1 


M3.8 


142.0 


3.1 


31 13.45 


0.006 


DD 


30 (2) 


-21 (2) 30 (3) 


-27 (3) 


296 273 


M2.4 


M3.4 


157.6 


6.0 


23 4.90 


0.015 


DD 


10(3) 


-40 (3) 5 (4) 


-40 (4) 


357 380 


M3.1 


M3.9 


163.6 


5.5 


22 2.26 


0.037 


DD 


-40 (2) 


-35 (2) -40 (2) 


-31 (2) 


242 270 


M2.2 


M3.1 


8.7 


3.9 


27 112.47 


0.005 


DD 


-48 (4) 


-5 (4) -41 (4) 


-3 (4) 


417 476 


M3.3 


M3.5 


174.1 


7.5 


58 1.59 


0.033 


DD 


37 (3) 


34 (3) 42 (4) 


31(4) 


429 445 


M1.5 


M2.3 


27.6 


5.6 


16 27.74 


0.034 


DD 


-9 (2) 


-42 (2) -7 (3) 


-41 (3) 


344 379 


M1.7 


M2.7 


11.5 


2.3 


34 74.69 


0.017 


DD 


36 (3) 


-16 (3) 37 (3) 


-21 (3) 


150 161 


M4.2 


M4.2 


17.0 


4.8 


11 21.02 


0.004 


DD 



Note. — The first 10 pairs are listed here; the full version of the table is available online. 
" The identifiers were generated using the standard Jhhmm±dd format using coordinates of the primary star and are prefaced with the string 'SLW. 

All magnitudes are psfmag and have not been corrected for extinction. Their errors are listed in parenthesis. Note that we use extinction-coiTected 
magnitudes in our analysis. 

The distances wer e calcu lated using photometric parallax relations and have 1 a errors of ~14%. The absolute distances to subdwarfs (SDs) are 
overestimated (see i|2.4. U . 

The spectral types were infeo'ed from the r — z colors jWest et al.l2008lJCovev et al.ll2007D and are correct to ±1 subtype. 
Class denotes the various types of pairs in SLoWPoKES. See Table |4] 

CPM triples are identified in our search (TableH]). Moreover, 
the current SLoWPoKES catalog is likely to contain quadru- 
ple systems in which the two components of the identified 
wide binary are themselves in fact spatially unresolved bina- 
ries with near-equal mass components. We have initiated an 
adaptive optics program to identify such higher-order systems 
in the SLoWPoKES catalog. 

4.L Kinematic Populations 

iLuvtenI (Il922h devised the reduced proper motion (RPM) 
diagram to be used in the stead of the H-R diagram when dis- 
tances to the objects are not available, as is the case in large 
imaging surveys. The RPM of an object is defined as 

H = m + 5 logM + 5 = Af + 5 legist - 3.25 (14) 

where ft is the heliocentric tangential velocity in km s^'^ 
given by vt = 4.74 i^i d and i^i is the proper motion in arcsec- 
onds yr^'^. Just as in a H-R diagram, the RPM diagram effec- 
tively segregates the v arious luminosity or kinematic classes 
interpretea by our algoritnm as navmg a aiscrepant pnoto- from each other (e.g. IChaname & Gouldll2004t iHarris et alJ 
metric distance from its wide tertiary companion. The avail- 120061: lOpine & Bongiornol2007l : ISIJ08h . 
able SDSS photometry and astrometry shows evidence of a WDs, in addition to being relatively very blue, are less lu- 

substantial number of such multiple systems, and we plan minous than either the DDs or the SDs; hence, the observed 
to make these the subject of a future study. Already, four WDs tend to be nearby disk WDs with high tangential ve- 



Table 4 

The SLoWPoKES binaries 



Class 


Type 


Number 


DD 


disk dwarf 


1245 


SD 


subdwarf 


70 


WD 


white dwarf-disk dwarf 


21 


T 


triple 


4 



CPM higher-order multiple systems in SDSS is very chal- 
lenging due to the lack of reliable proper motions in the 
SDSS/USNO-B matched catalog at Ad < 7" and r > 20. 
Unless all components are widely separated and are all bright, 
they will be rejected in our search. In addition, we are cur- 
rently rejecting hierarchical triples consisting of a close pair 
that is unresolved in SDSS and a wide, CPM tertiary. If the 
mass ratio of the close, unresolved pair is near unity, it will 
appear as an over-luminous single star that will then be mis- 
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Figure 8. 50" x 50"gri composite images of CPM pairs found in tlie SLoWPoKES survey. Pictured are liigli-mass ratio pairs (top row), identical twins (middle 
row), white dwarf-disk dwarf pairs (bottom row, left), and halo subdwarf pairs (bottom row, right). Spectral types, based on their r — z colors are shown. Overall, 
1 342 wide, low-mass binaries were identified. 

locities. Specifically, spectroscopic follow-up has shown that 
the disk WDs with vt 25-150 km s"^ (Figure [Tol black 
dashed lines) can be effectively identified from the 5-band 
RPM dia gram, when comp lemented by photometric parallax 
relations jKilic et all 12006). However, we note that spectro- 
scopic follow-up is needed to confirm that the identified ob- 
jects are actually WDs. The available SDSS spectra confirm 
that 9 of the 21 WD primaries, identified in SLoWPoKES, are 
indeed WDs. As the WDs were identified from the g-band, 
they scatter toward the SD and DD locus in the r-band RPM 
diagram. 

Subdwarfs are low-metallicity halo counterparts of the 
MS dwarfs found in the Galactic disk. Hence, they have 
bluer colors at a given absolute magnitu de (however, the 
g — r colors for M su bdwarfs are redder; iWest et aT] I2004t 
iLepine & Scholzir2008l) and have higher velocity dispersions. 
As a result, the subdwarfs lie below the disk dwarfs in the 
RPM diagranQ- To segregate the SDs from the DDs, we used 
the DD photometric parallax relations complemented by the 
mean tangential velocity of halo stars, vt = 180 km s^^. This 
relation is shown in blue dashed lines in the r-band RPM dia- 
gram in Figure [To] As the mean halo velocity was used, SDs 
can scatter above the line; but DDs would not be expected to 
be below the blue line. For comparison, vt ~ 25 km s^^, the 
mean tangential velocity of disk stars, is also shown. Note 
that DDs in SLoWPoKES have tangential velocities larger 
than the mean velocity of the disk, which is expected as we 
rejected all stars with i^l < 40 mas yr^^. 




AG (") 

Figure 9. The projected linear separation of SLoWPoKES pairs aie plotted 
as a function of angular separation, with the iso-distance lines plotted across 
the grid. DD (black), SD (red), and WD-DD pairs (blue) are shown. Sys- 
tematically lai'ger distances for the SDs is clearly seen, which is expected as 
photometric DD photometric parallax relations were used to calculate the SD 
distances as well. Hence, SDs are excluded in all analysis involving their 
distances. 



*^ ISIJ08I showed that the RPM diagram becomes degenerate at Z > 2-3 
kpc due to the decrease in rotational velocity. This does not affect our sample. 



14 



Dhital et al. 



12 
14 
16 

18F 
20 
22 




s_ , -v,.25 
' ^ v,=150 




-1 01234 01234 
g-i r-z 

Figure 10. SLoWPoKES pairs plotted in the Hg vs. g-i (left) and 
vs. r — 3 RPM diagrams; primaries (triangles) and secondaries (circles) of 
DD (black), SD (red), and WD-DD (blue) pairs are shown. The WDs are 
identified from the g-band RPM diagram, complemented by Mg(g — i) for 
WDs, as having 25 < ft < 150 km s~^ (black dashed lines) and are all 
expected to be part of the disk population. SDs are segregated from the DDs 
using the r-band RPM diagram, assuming SDs have vt > 180 km s^^ (blue 
dashed line) and the DD photometric parallax relations. 
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Figure 11. r vs. r — z Hess diagram and color-color plots for the SLoW- 
PoKES catalog; shown are the primary (triangles) and secondaries (circles) 
for DD (black), SD (red), and WD-DD pairs (blue). The subdwaifs are 
mostly of K spectral type and do not sho w redder q — r c olors at gi ven r — i 
ori — z colors, as seen in M subdwarfs jWest et ani2004i ; iLepine & Scholzl 



Figure [TOl shows the Hg vs. g — i and vs. 



RPM 



diagrams with both the primary (triangle) and the secondary 
component (circle) for the 21 WD-DD (blue), 70 SD-SD 
(red), and 1245 DD-DD (black) pairs that were identified in 
SLoWPoKES. The properties of the DDs, SDs, and WDs in 
various color-magnitude and color-color planes are compared 
in Figure[TT] The identified SDs are either K or early-M spec- 
tral types, as our magnitude and color limits exclude the M 
subdwarf locus. The overestimation in the distances to SDs 
is clearly evident in Figure |9] as they are at systematically 
larger distances relative to the DDs. As a result, the calcu- 
lated physical separations for the SDs are also systematically 
larger At present a substantial number of subdwarf candi- 
dates are rejected from SLoWPoKES because the overesti- 
mated distances result in Ad > 100 pc. 

RPM diagrams have been used to confirm the binar ity of 
CPM pairs in the past (e.g. IChaname & Gouldl |2004|) . As 
components of a binary system most probably formed from 
the same material and at the same time, they should be mem- 
bers of the same luminosity class (with the obvious exception 
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Figure 12. A randomly selected sample of SLoWPoKES pairs plotted on 
the RPM diagram; components of DD (dotted), SD (solid), and WD-DD 
(dashed) pairs are connected. As binaries are expected to be part of the same 
kinematic populations, they are expected to lie parallel to the dwarf/subdwarf 
tracks with the obvious exception of WD-DD pairs or when the two compo- 
nents lie within the error bars of each other 

of WD-DD pairs) and the line joining the components should 
be parallel to the track in which the systems reside. In the 
case of WD-DD dwarf systems or systems that have separa- 
tions comparable to the error bars in the RPM diagram, the 
line need not be parallel. The r-band RPM diagram, in Fig- 
ure[T2] confirms that the SLoWPoKES systems are associated 
pairs. 

4.2. Separation 

Despite the large number of optical companions found at 
larger angular separations, the final distribution of the identi- 
fied pairs is mostly of pairs with small angular separations. 
This is shown in Figure [T3] where, after a narrow peak at 
/S.0 ^ 15", the distribution tapers off and rises gently after 
A6' > 70". To convert the angular separations into physical 
separations, we need to account for the projection effects of 
the binary orbit. As that information for each CPM pair is 
not available, we apply the statistical correction between pro- 
jected separation (s) and true separation (a) determined by 
IFM92I from Monte Carlo simulations over a full suite of bi- 
nary parameters. They found that 



1.26 s = 1.26 Md 



(15) 



where the calculated a is the physical separation including 
corrections for both inclination angle and eccentricity of the 
binary orbit and is the actual semi-major axis. We emphasize, 
however, that these a values are valid only for ensemble com- 
parisons and should not be taken as an accurate measure of 
a for any individual system. In addition, the above equation 
implies that, at the extrema of angular separations probed, we 
are biased towards systems that are favorably oriented, either 
because of sinz projection effects or eccentricity effects that 
lead to a changing physical separation as a function of orbital 
position. For example, at A6' ~ 7", pairs at their maximal 
apparent separation are more likely to be identified while we 
are biased towards pairs with smaller apparent separations at 
A6I w 180". 

The distribution of a is plotted for the SLoWPoKES sam- 
ple in Figure [13] (right). The semi-major axes for the iden- 
tified pairs range from - 10^ - 10^ AU (~0.005-0.5 pc), 
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Figure 13. The projected angular se paiation (left) and the semi-major axes 
(right; corrected using a = 1.26 6* d lFM92D for the CPM doubles identified 
in SLoWPoKES. The upper x-axis in the right panel shows the semi-major 
axis in parsecs; the widest CPM pairs known in the literature (see Table|5) ai'e 
shown as pluses at their semi-major axis values at the top of the panel. SDs 
were excluded from the right panel. The distribution of physical separations 
is clearly bimodal with a break at ~0. 1 pc. 

with sharp cutoffs at both ends of the distribution, and with 
a clear bimodal structure in between. The cutoff at small 
separations is observational, resulting from our bias against 
angular separations A6 < 7". Indeed, the range of physi- 
cal separations probed by SLoWPoKES is at the tail-end of 
the log-n ormal di strib ution ( (a) = 30 AU, criog a ~ 1 -5) pro- 
posed by IDM91I and IFM921 The cutoff we observe at the 
other end of the distribution, a ^ 10^ AU, is likely physical. 
With the mean radii o f prestellar cores observed to be around 
0.35 pc or AU dBenson & Mverll ll989l: ICle mens et al.l 
IT99TFlJessop & Ward-Thompsonll2000h . SLoWPoKES repre- 
sents some of the widest binaries that can be formed. How- 
ever, other studies have found binary systems at similar sep- 
arations; plotted in pluses are the wide CPM pairs compiled 
from the Uterature and listed in Table |5] 



Between the cutoffs at ~ 10^ AU and at ^ 10^ AU, we ob- 
serve a distinct bimodality in the distribution of physical sepa- 
rations at a lO"*-^ AU (-^0.1 pc) that has no correlation with 
the distance to the observed system (see §|5]below). We have 
high confidence that this bimodality is not due to some sort of 
bias in our sample. As discussed above, the most important 
observational bias affecting the distribution is the bias against 
pairs with a < 10^ AU, because we are not sensitive to sys- 
tems with A9 < 7" nor to very nearby systems. In addition, 
given the care with which we have eliminated false positives 
in the sample, we have high confidence that the bimodal struc- 
ture is not due to a large contamination of very wide chance 
pairs. Instead, it is likely that this bimodality reveals two dis- 
tinct populations of wide binaries in SLoWPoKES, possibly 
representing systems that form and/or dissipate through dif- 
fering mechanisms. We discuss the bimodality in the context 
of models of binary formation and kinematic evolution in §|5] 

4.3. Mass Distribution 

Figure [14] (left panel) shows the r — z color distributions 
of the primary and secondary components of SLoWPoKES 
pairs, where the primary is defined as the component with the 
bluer r — z color (and, thus, presumably more massive). As 
the SLoWPoKES sample is dominated by dwarfs, the r — z 
color distribution should correspond directly to the mass dis- 
tribution. Both the primary and secondary distributions show 
a peak at the early-m id M spectral types (as inferred from 
their colors; § 12.4.2b . probably due to, at least in part, the 
input sample being mostly M0-M4 dwarfs. However, apart 
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Figure 14. Left: The r — z color for the primaries (solid) and secondaries 
(dashed). The upper x-axis, in the left panel, shows the spectral types inferred 
from their r ~ z colors £ovev et al. 2007; West et al. 2008). Both distribu- 
tions exhibit bimodality and are. similar to each other Right: The mass ratio, 
q = A/2 /Ml for SLoWPoKES pairs (soHd histogram) indicates that most 
are in equal- or near-equal-mass systems. Within the observational limits 
of SDSS, we should have been able to see much more extreme mass ratio 
systems (dashed histogram). The lack of such systems indicates that the ob- 
served distribution is not dominated by observational biases and is probably 
real. 

from the saturation at r < 14, there is no bias against find- 
ing higher-mass companions. Hence, it is notable that more 
than half of the primaries have inferred spectral types of MO 
or later, with the distribution of the secondaries even more 
strongly skewed to later spectral types (by definition). 

As discussed above, selection biases in SLoWPoKES are 
in general a function of distance because of the large range 
of distances probed (see Sec. 14. 11 Figures[3]and|9l). Thus, in 
Figure [15] we show the box-and-whisker plot of the distribu- 
tion of the r — z colors of the primary (blue) and secondary 
(red) components as a function of distance in 100 pc bins. 
The bar inside the box is the median of the distribution; the 
boxes show the inter-quartile range (defined as the range be- 
tween 25"^ and 75**^ percentiles); the whiskers extend to ei- 
ther the 1.5 times the inter-quartile range or the maximum or 
the minimum value of the data, whichever is larger; and the 
open circles show the outliers of the distribution. The black 
dashed lines show the bright and faint limits of our sample 
(14 ^ T < 20). As the figure indicates the primary and 
secondary distributions are bluer at increasing distances, as 
would be expected due to the bright and faint limits; but the 
secondary distribution does not change as compared to the 
primary distribution as a function of distance. However, this 
is likely to be a strict function of the faintness limit of our cat- 
alog: the stellar mass function peaks around M4 spectral type 
jBlOl) while we cannot see M4 or later stars beyond ^400 pc. 
Hence, with our current sample, we cannot discern whether 
the tendency for the secondary distribution to follow the pri- 
mary distribution is a tendency toward g ^ 1 or a tendency 
for the secondary distribution to be drawn from the field mass 
function. 

We can also examine the distribution of the secondary 
masses relative to their primaries, i.e. the mass ratio distri- 
bution, which is an important parameter in the study of bi- 
nary star formation and evolution. As shown in Figure [14] 
(right panel, solid histogram), the distribution of mass ratios, 
q = M2/M1, is strongly skewed toward equal-mass pairs: 
20.9%, 58.5%, 85.5% of pairs have masses within 10%, 30%, 
and 50%, respectively, of each other. To determine whether 
this is strictly due to the magnitude limits in SLoWPoKES, we 
calculated the mass ratios for hypothetical pairs with the ob- 
served primary and the faintest observable secondary. The re- 
sulting distribution (dashed histogram) is considerably differ- 
ent from what is observed and shows that we could have iden- 
tified pairs with much lower mass ratios, within the r < 20 
faintness limit. The same result is obtained when we pair 
the observed secondaries with the brightest possible primary. 
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Table 5 

A sub-sample of previously known wide binaries with the projected separation, s, > 10'* AU 



'0.05 pc) 





ID'' 


A0 
(") 


s 

(AU) 


Spectral Type'' 


BE" 
(10*'' ergs) 


References 


Primary 


Secondary 


Primary 


Secondary 


HIP 38602 


LSPM J0753+5845 


109.90 


10175 


G8.1 


M5.1 


132.18 


1 


HIP 25278 


HIP 25220 


707.10 


10249 


F1.8 


K5.8 


307.44 


1 


HIP 52469 


HIP 52498 


288.00 


10285 


A8.9 


F9.3 


379.74 


1 


HIP 86036 


HIP 86087 


737.90 


10400 


F8.4 


M2.6 


184.33 


1 


HIP 58939 


LSPM J 1205+ 1931 


117.20 


10558 


K0.5 


M4.7 


107.44 


1 


HIP 5 1669 


LSPMJ1033+4722 


164.30 


10668 


K7.3 


M5.7 


52.10 


1 


HIP 81023 


LSPMJ1633+0311N 


252.00 


10723 


K3.2 


M3.8 


81.51 


1 


HIP 50802 


LSPMJ1022+1204 


311.90 


11139 


K5.9 


M3.I 


56.44 


1 


HIP 78128 


LSPM J 1557+0745 


144.60 


11385 


G7.1 


M5.0 


123.01 


1 


HIP 116106 


2MASS J2331-04AB 


451.00 


11900 


F8 




161.10 


2 


a Cen AB 


NLTT 37460 


9302.00 


12000 


G2+K1 


M6 


430.21 


3 



N 




WBF = 



References. — (l)'LeDine & Bongiornd '20071); (2) ICaballerd '20071); (3) ICaballerd '200S»); f4) 'Caballerd '20101); (5) 
[Bahcall & Soneira (1981); (6) Latham et al. (1984); (7) Makarov et al. ( 2008); (8) Zapatero Osorio & Martin (20 0ll; (9 ) 
iFaheitv etalii201Ui) ; flO) .Chaname & Gould (.2004.) ; (ll) iOuiim et alji2009t) ; fl2) iPoveda et alil2009i) ; (13) i Allen et aimOOOl) . 

Note. — The first 10 pairs are listed here; the full version of the table is available online. 
^ We have tried to use HIP and NLTT identifiers, whenever they exist, for consistency. 

Spectral types are from the referenced papers, SIMBAD, or inferred from their V — J colors usin dKenvon & HartmanBlll99S . 

Binding energies are calculated using estimated masses as a function of spectral type IKraus & Hillenbranj|2007l) . When 
spectral type for the secondary was not available, it was assumed to be an equal-mass binary. 

WBF rises from ^0.23% at the bluest r — z color (~K7) to 
^0.57% atr — z = 1.6 ('^M2), where it plateaus. This trend 
is probably due at least in part to our better sensitivity to com- 
panions around nearby early-mid M dwarfs as compared to 
the more distant mid-K dwarfs or to the much fainter mid- 
late M dwarfs. To get a first-order measure of the effects of 
the observational biases, we can look at the WBF in specific 
distance ranges where all stars of the given r ~ z colors can 
be seen and the biases are similar for all colors, assuming the 
range of observable magnitude are r ~ 14-20. In particular, 
it would be useful to look at r — z = 0.7-1.5 where the WBF 
increases and r — z = 1 .5-2.5 where the WBF plateaus in Fig- 
ure[T6] Figure [it] shows the WBF in the two regimes for d = 
247-1209 pc and 76-283 pc, respectively, where all compan- 
ions of the given color range are expected to be detected at 
those distances. In the restricted ranges, the WBF is generally 
higher for a given color than in Figure [161 ranging between 
0.44-1.1%. More importantly, neither panel reproduces the 
trend in WBF seen for the same color range in Figure [16] in- 
dicating that observational biases and incompleteness play a 
significant role in the observed WBF. As these two samples 
are pulled from two different distance ranges, some of the dif- 
ferences in the value of WBF as well as the observed trend in 
WBF as a function of i — z might be due to Galactic structure 
or age of the sampled pairs. Hence, even the maximum ob- 
served WBF of 1.1% in SLoWPoKES is likely a lower limit 
on the true WBF. 

SLoWPoKES should be useful for studying the WBF as a 
function of Galactic height (Z), which can be taken as a proxy 
for age. Figure [TSlshows the WBF vs. Z,mr — z color bins. 
Again, to keep the observational biases similar across color 
bins, we selected distances ranges such that CPM compan- 
ions with r — z colors within ± 0.3 could be seen in SDSS; 
for example, in the first bin (r — z = 0.7-1.0) all companions 
with r — z = 0.4-1.3 are expected to be seen. Note that this 
approach is biased towards equal-mass pairs, with high mass- 
ratio pairs never counted. Consequently, the WBF is lower 
than in the Figures [16] and [17] with maximum at ~0.35% for 
r — z = 1.3-1.9; the WBF peaks at around the same color 
range in all three figures. However, our method ensures that 
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Figure 15. The color distribution for the primary (blue) and secondary (red) 
components of SLoWPoKES pairs as a function of distance, in 100 pc bins 
shown in a box-and-whisker plot; the number of pairs in each bin is printed 
along the top. The bar inside the boxes show the median of the distribu- 
tion; the boxes show the inter-quartile range; the whiskers extend to either 
1.5 times the inter-quartile range or the maximum or the minimum value of 
the data, whichever is larger; and the open circles show the outliers. The 
dashed lines show the magnitude limits of this catalog (14 < r < 20). The 
secondary distribution tracks the primary distribution at all distances. 

Hence, we conclude that the observed distribution peaked to- 
ward equal-masses among wide, low-mass stars is real and 
not a result of observational biases. 

4.4. Wide Binary Frequency 

To measure the frequency of wide binaries among low-mass 
stars, we defined the wide binary frequency (WBF) as 

mimbcr of CPM pairs found 



number of stars searched around 



(16) 



With 1342 CPM doubles among the 577459 stars that we 
searched around, the raw WBF is 0.23%. However, given 
the observational biases and our restrictive selection criteria, 
0.23% is a lower limit. Figure [TSlshows the WBF distribu- 
tion as a function of r — z color with the primary (solid his- 
togram), secondary (dashed histogram), and total (solid dots 
with binomial error bars) WBF plotted; the total number of 
pairs found in each bin are also shown along the top. The 
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Figure 16. The wide binary fraction (WBF) of low-mass stars that are the 
primary (solid histogram) and secondary (dashed histogram) companions in 
wide pairs shown as a function of their r — z colors, with spectral types 
shown along the top x-axis. The total WBF of low-mass stars in wide pairs 
, shown in black circles along with the binomial errors, shows a peak of ~ 
0.57% for r — z = i.6 (~M2). The low WBF is expected due to the severe 
observational biases, as well as a reflection of our restrictive binary selection 
algorithm. The total number of CPM pairs in each bin are printed along the 
top. 
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Figure 17. Same as Figure [T6l but to explore the effect of biases in the WBF, 
we compare two color ranges: r — z = .7- 1.5 where the WBF rises and 
r — z = 1.5-2.6 where it plateaus in Figure [T6l To make the selection effects 
similar throughout the selected color range, we restricted the distance of the 
stars such that all secondaries of the given color range can be seen throughout 
the distance range. The WBF seen here (i) is higher than in Figure 1 1 6[ for a 
given color bin an d (i i) as a function of r — 2 does not reproduce the same 
trend as in Figure [16] Since the two samples were chosen from two differ- 
ent distance ranges, they might also reflect differences due to the age of the 
sampled pairs and/or the Galactic structure at that position. These differences 
indicate significant observational biases and incompleteness affect our WBF 
determination. 

we are sensitive to all similar mass pairs across all of the 
distance bins. Figure [18] suggests that the WBF decreases 
with increasing Z. As this trend appears for both primary and 
secondary components for almost all spectral types (colors) 
probed, it is likely not an artifact of observational biases and 
is strong evidence for the time evolution of the WBF. Wide bi- 
naries are expected to be perturbed by inhomogeneities in the 
Galactic potential, giant molecular clouds, and other stars and, 
as a result, to dissipate over time. As pairs at larger Galactic 
heights are older as an ensemble, it is expected that a larger 
fraction of pairs at larger Z, which are older, have dissipated. 



5. DISCUSSION 

As discussed above, ISIJ08I recently conducted a search 
of CPM pairs in the SDSS DR and found -22000 wide 
pairs, of all spectral types, with the probability of being a 
real binary of —67% at angular separations of 5-30". This 
search was done for all stars in the DR6 photometry with 
fj, > 15 mas yr~^ and match ed both components of proper mo- 
tions within 5 mas yr~^. As lSIJOSl note. the sample is not very 
efficient for follow-up studies due to the hkely large number 
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Figure 18. WBF, with binomial errors, decreases with Galactic height at all 
r — 2 color bins for the primaries (circle), secondaries (triangle), and the total 
(solid histogram). In order to keep biases similar in a color bin, we selected 
a distance range such that all stars with r — z colors within ±0.3 dex of 
the bin could be seen throughout the distance range. This decrease in WBF 
with increasing Z suggests dynamical destruction of (older) pairs at larger 
Galactic heights. 

of chance alignments (i.e. false positives). Our analysis in this 
paper indicates that their low proper-motion cutoff is chief 
culprit for the contamin ation; w ith a higher proper-motion 
cutoff similar to ours, the iSIJOS] sam ple could in principle be 
sifted of many false positives. IsiJOSi used a statistical subtrac- 
tion to correct the ensemble for the chance alignments, which 
suffices for the purpose of their work. While SLoWPoKES is 
intended to serve a very different purpose — namely to provide 
a "pure" catalog of high-fidelity systems suitable for targeted 
follow-up studies — it is useful to compare so me of o ur tenta- 
tive findings and interpretations with those of lSIJOSl based on 
their more complete sample. 

For example, the u ppe r limit of the WBF in SLoWPoKES 
(—1.1%; e.g . Figur es [ITIi is consistent with the 0.9% at Z = 
500 pc from lSIJOSl This suggests that, for specific color and 
distance ranges, the selection biases in SLoWPoKES may 
not significantly affect the ability of this sample to charac- 
terize the WBF among late - type s tars at the widest separa- 
tions. iKraus & Hillenbrand! (12009^ similarly find a WBF of 
at most a few percent for masses 0.012-2.5 Mq at separa- 
tions of 500-5000 AU in their study of the Taurus and Up- 
per Sco star-forming regions. Moreover, they find that the 
distribution of binary separations in their sample remains flat 
out to the limits of their survey (—5000 AU) for all but the 
lowest mass systems with Mtot< 0.3 Mq. At the same 
time, all of these studies fall well below the WBF (for s > 
1000 AU) of —9.5% determine d for solar-type Hipparcos 
stars (Opine & Bongiornoll2007[) and the 10% suggested by 
iLonghitano & BinggeUl (120 lOt) for G5 or later stars in the so- 
lar neighborhood from SDSS D R6. Similarly, integrating the 
Gaussian distribution of lDM911 suggests a WBF of —9.2% for 
1000 AU < a < 0.5 pc and -0.64% for a > 0.5 pc for solar- 
type stars. E ven acc ounting for the lower multiplicity seen in 
M dwarfs bv lFM92l this is much higher than our result. Simi- 
larly, the decrease of WB F with increasing Z in SLoWPoKES 
(Figure [TSTi was noted bv lSIJOSl in their sample as well. 

S1J08 had noted that the companions of red stars, unlike the 
bluer stars, seemed to be drawn randomly from the field stellar 
mass function. We do not find evidence to support that find- 
ing, but our results are probably limited by the faintness limit. 
With respect t o the distribution of mass ratios (Figure [TSt , 
[Reid & GizisI (Il997h observed, as we do, a strong skew to- 
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wards equal-mass pairs w ith a peak at q >0.8 in their 8 pc 
sample. However. IFM92I found a flat dis tribution in their sam- 
ple of dMs within 20 pc. More recently jKraus & Hillenbrandl 
^200^^1) also found that a bias towards equal-masses in their 
study of 0.012-2.5 M© wide binaries (500-5000 AU sep- 
arations) in the Taurus and Upper-Sco star-forming regions. 
In contr ast, solar-type field stars exhib it a definitive peak at 
q ~ 0.3 (lDM9lLlHalbwachs et al.ll2003l) . With the field mass 
function peaking around '^M4 spectral type (IBIOI) . a peak at 
q ~0.3 for solar-type stars and alq^l for M dwarfs are both 
consistent with the secondary mass function being a subset of 
the field mass function. 

With some systems wider than ^ 1 pc, SLoWPoKES pro- 
vides the largest sample to date of the widest CPM doubles. 
However, these are not the only systems known at this separa- 
tion regime; Table |5] lists 84 systems with projected physical 
separations greater 0.05 pc (^lO** AU). This sample is biased 
towards pairs with at least one Hippa rcos star, with about 
half o f them from the CPM catalog of iLepine & Bongiornol 
(120071) . This is not a comprehensive list of all kno wn wide 
binaries; for example, the existing CPM catalogs of iLuytenI 
^Wf§. T988^. 'Bahcall & Soneiral (UMl). iHalbwachsl (Il986h. 
and Chaname & Gould (2004) likely contain at least a few 
hundred pairs in this separation range. Even among the 
brown dwarfs, where most known pairs have separations 
smaller than 15 AU, two very-wide systems— the 5 100 AU 
2MASS J0126555-502239 dArtigau et al."2007','2009^ and the 
6700 AU 2MASS J12583501H-40130 (Radigan et al. 2009)— 
have already been identified. All of these wide pairs and 
other low -binding energy VLM systems from the VLM binary 
archiveE3 are shown, when relevant, in Figures [T3]andl20lfor 
comparison with the SLoWPoKES distribution. Most of the 
CPM doubles in Table |5] were found in nearby, high proper 
motion catalogs. Due to the depth of SDSS imaging and the 
larger distances probed, SLoWPoKES significantly increases 
CPM pairs with large physical separations. For the same rea- 
son, the existence of pairs as wide as > 1 pc in SLoWPoKES 
(Figure [T3Tl should perhaps not be surprising. 

At the same time, the question of how such wide, weakly 
bound systems survive over time and how they form in the 
first place is a very interesting issue that the SLoWPoKES 
sample will be well suited to address. The initial distribution 
of separations of wide binaries is not static over time but is 
modified by interactions with other stars, molecular clouds, 
and variations in the Galactic potential. Over the lifetime of 
a binary, the small and dissipative, but numerous, encounters 
that it undergoes with other stars is much more efficient at dis- 
rupting the s ystem than single cata strophic encounters, which 
are very rare (IWeinberg et al.l[T987l) . Using the Fokker-Planck 
coefficients to describe the advection and diffusion of the or- 
bital binding energy due to such small encounters over time, 
[Weinberg et al. (1987|) calculated that the average lifetime of 
a binary is given by: 
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where and A/^ are the number density and average mass 
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Figure 19. The physical separation of SLoWPoKES pairs (black circles) and 
previously known wide CPM doubles (red pluses; Table[5) and VLM binaries 
(blue pluses) as a function of the tota l mass of the syste m, as infen'ed from 
thei r r — z colors. Previou s studies bv lReid et alj I200"T^) (green dashed line) 
an d IBurgasser eTal] f2003l) ( green dotted line) have also suggested empirical 
limits for stability. As neither of those describe the SLoWPoKES envelope, 
we have followed these previo us authors and fit a log -normal (solid green 
fine). Numerical simulations bv lWeinberg et alj 119871) suggest that the life- 
time of wide pairs is a function of age as well; the "isochrones" show the 
expected a vs. Mtot relationships for dissipation times of 1, 2, and 10 Gyrs 
(solid red lines), as defined by Eq. (Ts). While some very wide SLoWPoKES 
pairs are expected to dissipate on timescales of 1-2 Gyr, most should be sta- 
ble for > 10 Gyrs. All SDs are excluded from this figure. 

of the perturbers, V^ci is the relative velocity between the bi- 
nary system and the perturber, il/tot and a are the total mass 
and semi-major axis of the binary, and A is the Coulomb 
logarithm. Using the average Galactic disk mass density of 
0.11 Mq pc~'^, an average pertur ber mass of 0.7 Mp ,, Vtci = 
20 km s^^, and A = 1 (see also IClose et al.ll2007l) . we can 
rewrite the above equation as; 



1 mo -^tot 

a ~ 1.212 pc 



(18) 



for M in M0 and t^, in Gyrs. This describes, statistically, 
the widest binary that is surviving at a given age. For ex- 
ample, at 1 Gyr 1 M0 binaries as wide as 1.2 pc are likely 
to be bound, but by 10 Gyr, all systems wider than 0.12 pc 
will likely have been disrupted. Hence, the combination of 
Eq. (fTSl l and the distribution of binary separation at birth may 
describe the current binary population of the Galaxy. Fig- 
ure[T9]shows the physical separation vs. total mass for SLoW- 
PoKES pairs and the other known wide CPM pairs from the 
literature (see above); characteristic disruption timescales of 
1, 2, and 10 Gyr for Eq. (fTSl l are over-plotted. From this, 
most SLoWPoKES pairs can be expected not to dissipate be- 
fore 1-2 Gyr, with approximately half of the population stable 
enough to last longer than 10 Gyr. 

Interestingly, the widest SLoWPoKES systems appear to 
violate previously proposed empirical limits on maximum 
separations (dashed and dotted g reen lines in Figure [T9l 
iReid et al.l 12001 at iBurgasser etaITl2003h . Similarly, previ- 
ously proposed empirical limits based on binding energies, 
10^^ erg s for stellar and and 10"^^-^ ergs very low- mass 
regimes (IClose et al.l 120031 120071; IBurgasser et al.ll2007h . are 
also violated by SLoWPoKES pairs as well as other known 
CPM doubles. This is clearly evident in Figure|20l where we 
see that the most weakly bound SLoWPoKES systems have 
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binding energies of only ^ 10 ergs, comparable to the bind- 
ing energy of Neptune about the SurPl. iFahertv et al.l ( 120101) 
have also noted such transgressions in their sample of wide, 
very low-mass CPM doubles. 

It seems that previously proposed empirical limits were too 
restrictive as the widest and the VLM pairs that clearly vio- 
late them w ere not available a t the time. We can follow the 
approach of lReid et all (1200 lal) and lBurgasser et alj (120031) to 
fit a log-normal that forms an empirical envelope around the 
SLoWPoKES pairs (Figure [191 green solid lines). Although 
this new fit is not likely to define the absolute stability limit 
for wide pairs, the fit is notable for how well it describes the 
envelope. 

We note that Eq. (fTsT l as well as most previously proposed 
empirical limits have used the total mass of the system and do 
not take into account the mass ratio of the components. This 
creates a degeneracy as the binding energy is dependent on 
the product of masses. Instead, it may be more physical to 
consider the wide, low-mass pairs as a function of the product 
of the masses of their components — i.e., the binding energy 
(BE) of the system. Figure |20] (left panels) present the a and 
BE vs. Ml X A/2 for the SLoWPoKES pairs. That the disti-i- 
bution of a as a function of either sum or product of the com- 
ponent masses appears qualitatively similar is a manifestation 
of the fact that SLoWPoKES pairs tend to be equal mass (see 
Figure [T4]i: for high-g pairs there would be a more noticeable 
difference in these distributions. 

Perhaps, of particular interest to the issue of wide binary 
stability, formation, and evolution is our observation of a dis- 
tinctly bimodal distribution of binary separations, seen in Fig- 
ure 1^ which equates to a bimodal distribution of system 
binding energies as well. The break in the bimodal distri- 
bution at binding energy ~10^^'^ ergs is remarkable for how 
well it agrees with the previously proposed empirical limit in 
binding energy for stellar binaries to be dyna mically stable 
(IClose et al.ll2003l l2007t iBurgasser et"ani2007h . In addition, 
as can be seen in F i gure[l9l the 10 Gyr disruption "isochrone" 
of [Weinberg et al.l ( 1987b nicelv traces the observed break in 
the distribution over nearly a full decade of total system mass. 
The trend of the break in the distribution with mass is also im- 
portant because it makes it very unlikely that the bimodality is 
the result of observational biases in the SLoWPoKES sample. 
Similarly, there is no effect of the heliocentric distances or the 
Galactic heights (Figure |20l right panels) of the pairs on the 
physical separation or binding energy. As discussed earlier, 
the principal observational biases in SLoWPoKES affect only 
(i) the smallest separations most severely because of the strict 
cut-off at 7" and (ii) to a lesser extent the largest separations 
as we limit our search to 180". 

Thus, the bimodal distribution of binary separations sug- 
gests two populations of binaries, perhaps representing (i) 
systems of stars that formed with sufficient binding energy to 
survive for the age of the Galaxy and (ii) relatively young sys- 
tems that formed within the past 1-2 Gyr but that will likely 
not survive much longer. Figure |2T] shows that the dispersion 
in the tangential velocity is indeed smaller for the "young" 
population, relative to the "old" population that has had more 
time to get kinematically heatedH Quantitatively, the median 

SLoWPoKES systems are much more hkely than the Sun-Neptune sys- 
tem to be disrupted by passing stars due to the binaries' larger cross-section 
of interaction. 

" One might expect lower mass stars to move with higher velocities rel- 
ative to higher mass stars if equipartition of kinetic energy holds. However, 



absolute deviations for the two populations are 46.1 km s~ 
and 32.7 km s~^. Curiously, the separations of the binaries in 
the second group (up to 1 pc; Figure [13] right panel) is larger 
than the typical sizes of prestellar cores (~0.35 pc), suggest- 
ing that these systems may not have formed via the "standard" 
collapse and fragmentation of individual cloud cores. 

Recent N-body simulations predict bimodal distributions 
of wide binary separations and may provide important the- 
oretical insight for understanding the formation and evolu- 
tion of wide pairs. Although the specific simulations pre- 
dict bimodality on different scales than seen in our observa- 
tions, some of the resulting predicted physical properties of 
pairs are similar to those fo und in SLoWPoKES. For example, 
iKouwenhoven et aTj ( 120101) find that the distribution of binary 
semi-major axes will in general be bimodal on the size scale 
of star-forming clusters (i.e. '^0.1 pc), with a group of tight 
systems formed primarily via dynamical interactions near the 
cluster center and a second group of loosely bound systems 
formed through random pairing during the slow dissolution 
of the cluster In their simulations, up to ^ 30% of the re- 
sulting binaries will be in the latter "dissolution peak", with 
separations of ~ 0.1 pc. They predict that the mass ratios 
of these systems will refl ect random pairing. Alternatively, 
iJiang & Tremainel (l2009h have found that a bimodal distri- 
bution arises due to the slow (as opposed to instantaneous), 
diffusive disruption of systems on timescales of several Gyrs, 
and that the details of the resulting bimodality is mostly inde- 
pendent of the initial distribution of separations arising from 
the formation process. Their simulations predict that the bi- 
modal break in the distribution occurs at a few Jacobi radii, 
which corresponds to a separation of ~ 2 pc for a total sys- 
tem of 1 Mq . Perhaps the diffu sive dissipation mechanism is 
more efficient than predicted bv lJiang & Tremai"n3 (120091) . or 
perhaps the "tidal-tail peak" of dissolving systems has not yet 
been observed. In the future, the SLoWPoKES catalog can 
be extended to even larger physical separations to probe this 
possibility. 

6. CONCLUSIONS 

We have created the SLoWPoKES catalog, comprising 
1342 CPM binary pairs, identified through statistical match- 
ing of angular separation, photometric distances, and proper 
motion components. We have sifted the sample of chance 
alignments using a Galactic model based on empirical obser- 
vations of stellar spatial and kinematic distributions. With the 
objective that each pair can be confidently used to investigate 
various science questions regarding low-mass stars, we have 
adopted a very restrictive set of selection criteria. This ap- 
proach clearly underestimates the number of binary systems. 
Moreover, the sample includes several biases, the most impor- 
tant of which are a lack of systems with physical separations 
smaller than ~ 1000 AU (arising from a strict bias against an- 
gular separations smaller than 7"), and the exclusion of cer- 
tain types of higher-order multiples (e.g. triples) due to the 
strict photometric distance matching. However, as a conse- 
quence the catalog should contain very few false positives, 
making follow-up studies efficient. 

We built a Monte Carlo-based six-dimensional Galactic 
model that is able to replicate the positional and kinematic 

observations have shown that the equipartition of kinetic energy does not hold 
and that kinematic heating is purely a function of age for field stars. Unlike 
in clusters, individual stars can be thought of as small, point masses when 
compared to the much larger Galactic potential or the giant molecular clouds 
that cause dynamical heating dEoc hanski et al. 2007; Schmidt et al. 2dTcl) . 
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Figure 20. The physical separations and binding energies of SLoWPoKES pairs as a function of the product of their masses (left) or the Galactic height (right), 
with the histograms of the distributions shown along the top and the right. For comparison, other CPM pairs (red pluse s; Ta blets) and VLM binaries (blue pluses; 
VLM binary archive) from the literature are shown. The bimodality in physical separation (at ~0.1 pc) seen in Figures fTsl and figl and along the top panels in this 
figure is evident in binding energy as well (at r^W^^ '^ ergs); these are marked with blue dashed lines in the histograms. As the binding energy is not dependent 
on Z, any selection effects or biases that depend directly or indirectly on Z are not likely to be the cause of the bimodality. Previously suggested minimum 
binding energy for stellar (lO'*^'^ ergs) and substellar (10'*" ergs) binaries IClose et al.l2003l) are violated by most SLoWPoKES as well as other CPM and VLM 
pairs. All SDs are excluded from this figure. 

is how difficult it is to find two physically unassociated stars 
close together in space: along a typical SDSS LOS, there are 
expected to be only 0.52 chance alignments within 15" and 
a minuscule 0.03 chance alignments if the volume is consid- 
ered. The additional matching of proper motions gives each 
of the accepted SLoWPoKES binaries a very low probability 
of being a false positive. 

Due to their intrinsic faintness and the resulting small num- 
bers, binarity studies of low-mass stars have been limited in 
scope. However, with the advent of large-scale deep sur- 
veys, detailed and statistically significant studies of M (and 
L) dwarfs are being done. SLoWPoKES is now the largest 
sample to date of very wide, low-mass binaries. In particular, 
SLoWPoKES provides a large sample of systems with phys- 
ical separations up to a 1 pc that will be useful for putting 
firm constraints on the maximum size of physically associ- 
ated systems. How the widest of these systems form, and 
how long they survive, is in particular an interesting question 
that SLoWPoKES is well suited to address. While numerical 
calculations suggest that approximately half o f SLoWPoKES 
system s can remain bound for at least 10 Gyr ([Weinberg et alj 
119871) . previously proposed empirical limits are violated by 
many SLoWPoKES systems. 
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Figure 21. The tangential velocity (tit = 4.74 /i d) distribution for the 
low-binding energy population (dashed histogram) and high-binding energy 
distribution (solid histogram), with the segregation at binding energy of 
^q41.3 ergs that was noticed in Figure|20] A smaller dispersion in the tangen- 
tial velocity supports our suggestion that the low-binding energy distribution 
is of a younger population. 

properties of the stars in the MiUcy Way. In its current in- 
carnation, we used it to calculate the the number of stars 
within a certain spatial volume in the Galaxy and the likeli- 
hood that those stars have common kinematics (proper mo- 
tions) by chance. One of the things this model underscores 
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Indeed, the distribution of SLoWPoKES binary separations 
is distinctly bimodal, suggesting the presence of (i) a popula- 
tion of tightly bound systems formed with sufficient binding 
energy to remain intact for the age of the Galaxy and (ii) a 
population of weakly bound systems that recently formed and 
that are unlike ly to survive past 1-2 Gyr Recent N-body sim- 
ulatio ns (e.g. iKouwenhoven et all 120 lOt iJiang & Tremaind 
120091) in fact predict a bimodal distribution of binary sepa- 
rations on the scales probed by the widest SLoWPoKES sys- 
tems. 

We observed a wide binary frequency of ^1.1%, which is 
likely to be a minimum given the natu re of ou r sample. While 
this is consistent with the results from lSIJOSi who found 0.9% 
of stars at Z = 500 pc had wide companions., it is signifi- 
cantly lower than ~9.5% of nearby solar-type Hipparcos stars 
having wide companions (Lepine & Bongiorno 2007). While 
the incompleteness involved in the searching for companions 
at la rge dis tances probably causes some of this, both this study 
and lSIJOSi saw a decrease in binary fraction as a function of 
Galactic height, evidence of dynamical destruction of older 
systems. Hence, the wide binary fraction around our initial 
sample of low-mass stars might actually be significantly lower 
than the Hipparcos stars. 

Besides the importance of SLoWPoKES for constraining 
models of formation and evolution of binary stars, SLoW- 
PoKES systems are coeval laboratories — sharing an identi- 
cal formation and evolutionary history without affecting each 
other — making them ideal for measuring and calibrating em- 
pirical relationships between rotation, activity, metallicity, 
age, etc. We have started programs t o test and calibrat e the 
age-activity relationship measured bv IWest et al.l ( 2008 ) and 
to explore whether gyrochronology (lBarnesll2003l 2007 ) can 
be applied in the fully-convective regime. As coeval labora- 
tories allow for the removal of one or more of the three fun- 
damental parameters (mass, age, and metallicity), much more 
science can be done with a large sample of such systems. 

Future astrometric missions, such as the Space Interfer- 
ometry Mission (SIM), should provide exquisite astrome- 
try, perhaps enabling us to trace the orbits of some of the 
SLoWPoKES systems. While tracing orbits with periods 
> years sounds ambiti ous, with SIM's m icroarcsec- 

ond level (or better) astrometry (lUnwin et al.ll2008l) combined 
with SDSS, DSS, and/or other epochs, it is not unrealistic. 
Similarly, the "identical" twins in SLoWPoKES would be 
ideal sites to probe for the presence and differences in the 
formation mechanism of planets. As each identical twin in 
a CPM double provides similar environment for the forma- 
tion and evolution of planets, these systems can be ideal sites 
to study planetary statistics. Due to their large separations the 
stars are not expected to influence each others evolution but 
have similar mass, age, and metallicity, as we noted earlier in 

The SLoWPoKES catalog, as the name suggests, only con- 
tains systems for which kinematic information is available. 
We can, however, use the results from the Galactic model to 
identify pairs at the small separations {/S.9 < 7"), albeit with 
a larger uncertainty, for which no kinematic information is 
available. Similarly companions which are fainter than r = 
20 can also be identified as the SDSS photometry is complete 
to r = 22.5. The latter systems are likely to be skewed to- 
wards late-type dMs and unequal-mass pairs. A follow-up 
paper will study such systems and will add a large proportion 
of wide systems. 
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